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ABSTRACT 

This  paper  treats  the  excitation  of  a  transverse  magnetic  surface 
wave  on  a  dielectric  cylinder.   Specifically,  the  efficiency  with  which 
a  circular  filament  of  magnetic  current  excites  the  E   mode  on  an 
infinite,  dielectric  rod  is  determined.   The  E   mode  is  the  lowest 
order,  transverse  magnetic  mode  which  can  propagate  on  a  dielectric  rod. 
The  excitation  efficiency  of  a  source  is  defined  as  the  ratio  of  the 
power  converted  to  the  surface  wave  mode  to  the  total  power  which  is 
delivered  by  the  source. 

The  solution  for  the  field  generated  by  the  magnetic  current  fila- 
ment is  obtained  in  the  form  of  an  inverse  Fourier  transform.   The  inte- 
gral is  evaluated  by  considering  it  as  a  contour  integral  in  the  plane 
of  the  complex  propagation  constant.   The  far  zone  radiation  field  of 
the  source  is  obtained  from  the  asymptotic  evaluation  of  the  integral 
by  means  of  the  saddle  point  method  of  integration.   The  poles  of  the 
integrand  yield  the  TM  surface  wave  modes,  which  are  all  circularly 
symmetric.   By  restricting  the  range  of  k  b,  which  is  the  circumference 
of  the  dielectric  rod  in  wavelengths,  it  is  assured  that  only  the  E 
mode  is  launched  by  the  source.   Curves  are  presented  which  show  excita- 
tion efficiency  as  a  function  of  k  a,  the  circumferential  length  of  the 
filament.   It  is  found  that  a  filament  of  the  proper  diameter  will 
excite  the  E_  mode  with  an  efficiency  of  approximately  95  percent. 

In  the  laboratory,  a  very  narrow,  annular  slot  in  a  large  metal 
sheet  was  used  to  approximate  the  magnetic  current  filament.   The  excita- 
tion efficiency  of  the  slot  was  measured  using  Deschamps'  method  for 
calibrating  a  two-port  waveguide  junction.   The  measured  efficiencies  are 
in  close  agreement  with  the  results  predicted  by  the  theory. 
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1 .   INTRODUCTION 

A  number  of  papers  published  in  the  last  decade  have  treated  the 
mode  characteristics  of  surface  waves  on  various  types  of  open  wave- 
guides.  Of  particular  interest  has  been  the  utilization  of  such 
structures  as  surface  wave  antennas .   A  factor  of  prime  importance  in 
all  surface  guide  applications  is  the  efficient  excitation  of  the 
desired  mode  on  the  guide.   The  excitation  efficiency  of  a  source  is 
defined  as  the  ratio  of  the  power  converted  to  the  surface  wave  mode 
to  the  total  power  which  is  delivered  by  the  source. 

Most  of  the  published  articles  on  surface  wave  excitation  efficiency 
have  treated  sources  placed  near  infinite  planar  structures  such  as  the 
dielectric  coated  metal  sheet.   In  order  to  determine  the  excitation 
efficiency  of  a  source,  one  must  solve  the  source  form  of  Maxwell's 
equations  which  amounts  to  solving  an  inhomogeneous  wave  equation.   The 
usual  technique  of  solving  the  differential  equation  is  to  apply  the 
method  of  integral  transforms  which  yields  the  solution  for  the  field 
in  the  form  of  a  definite  integral.   The  integral  is  evaluated  by 
considering  it  as  a  contour  integral  in  the  plane  of  the  complex  propa- 
gation constant.   The  integrand  has  poles  and  branch  points  in  the  com- 
plex plane.   By  deforming  the  contour  of  integration  in  a  suitable  man- 
ner, the  integral  becomes  equal  to  the  sum  of  the  residues  at  the  poles 
plus  a  branch  cut  integral.   The  poles  correspond  to  the  surface  wave 
modes  and  the  amplitude  of  a  surface  mode  is  given  by  the  residue  of 
the  integrand  at  the  pole.   The  branch  cut  integration  yields  the 
radiation  field.   Alternatively,  the  radiation  field  is  obtained  from 
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the  asymptotic  evaluation  of  the  integral  by  means  of  a  saddle  point 
integration.   This  method  of  analysis  yields  the  amplitude  of  the 
radiation  field  and  of  the  surface  mode.   One  may  then  calculate  the 
radiated  power  and  surface  wave  power  and  determine  the  excitation 
efficiency  of  the  source. 

Several  authors  have  applied  this  method  to  particular  problems. 
Cullen   calculated  the  efficiency  of  an  infinitely  long  slot  above  a 
dielectric-coated  plane  conductor.   He  also  treated  the  slot  above  a 
corrugated  surface  and  the  case  of  a  chopped  surface  wave  source  such 
as  a  gently  flared  horn.   His  theory  predicted  efficiencies  greater 

than  90  percent  for  the  infinite  slot  above  the  corrugated  surface. 

2 
Whitmer  determined  the  efficiency  of  an  infinitely  long,  electric 

current  filament  placed  in  a  dielectric  sheet.   The  maximum  efficiency 

calculated  for  this  source  was  approximately  85  percent.   The  problems 

treated  by  Cullen  and  Whitmer  are  primarily  of  academic  interest,  since 

they  involve  infinite  source  distributions  and  doubly  infinite  planar 

structures.   As  yet,  the  efficiencies  predicted  in  Cullen' s  and  Whitmer3 s 

work  have  not  been  verified  experimentally. 

3,4 
Roberts    investigated  the  single  wire  transmission  line  excited 

from  a  flanged  coaxial  cable.   This  source  launches  the  transverse 

magnetic,  E   mode  on  the  wire.   He  evaluated  the  input  conductance  of 

the  coaxial  line  by  assuming  an  infinitesimal  gap  between  the  wire  and 

the  outer  conductor  of  the  coax.   He  defines  a  radiation  conductance  G 

and  a  surface  mode  conductance  Gn.   The  excitation  efficiency  of  the 

source  may  be  obtained  from  G  and  G_.   Roberts  measured  the  properties 

of  two  different  wires  in  the  laboratory.   Part  of  his  experimental 


3 
investigation  included  measuring  the  input  conductance  of  the  coaxial 
line  for  a  number  of  different  gap  radii.   The  measured  conductances 
were  plotted  on  a  graph  and  the  curve  was  extrapolated  to  obtain  the 
conductance  for  an  infinitesimal  gap,  since  it  was  impossible  to 
measure  such  a  source  experimentally.   The  extrapolated  value  agreed 
quite  closely  with  the  conductance  predicted  from  theory .   To  the 
author's  knowledge,  this  is  the  only  previous  experimental  corrobor- 
ation of  efficiency  as  predicted  from  the  pole,  branch  cut  analysis. 

Since  only  limited  experimental  verification  of  this  method  of 
analysis  has  appeared  in  the  literature,  this  research  was  directed 
toward  obtaining  more  complete  experimental  confirmation  for  a  practi- 
cal theoretical  problem.   The  surface  waveguide  which  was  chosen  for 
the  problem  is  the  dielectric  rod  since  it  is  a  surface  wave  structure 
which  has  practical  application.   The  source  is  a  circular  filament  of 
magnetic  current  placed  inside  the  dielectric  rod  concentric  with  the 
longitudinal  axis  of  the  rod.   The  magnetic  current  ring  was  selected 
since  it  is  a  source  which  has  a  simple  physical  realization.   One 
would  expect  the  ring  to  be  an  efficient  exciter  of  the  lowest  order, 
circularly  symmetric,  transverse  magnetic  mode  which  can  propagate  on 
a  dielectric  rod.   This  mode  is  known  as  the  E_   mode.   The  purpose  of 
this  investigation  is  to  determine  the  efficiency  with  which  the  magnetic 
current  ring  excites  the  E   mode  using  the  pole,  saddle  point  integration 
method,  and  then  to  measure  the  efficiency  experimentally.   It  will  be 

seen  later  that  excitation  efficiency  may  be  readily  measured  using 

5,6 
Deschamps     method  for  determining  the  scattering  matrix  coefficients 

of  a  two  port  waveguide  junction. 
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A  similar  problem  has  been  treated  by  C.  Jauquet ,   who  selected  a 

magnetic  current  ring  which  was  greater  in  diameter  than  the  dielectric 
rod.   In  his  paper,  Jauquet  obtains  an  integral  solution  for  the  field 
but  does  not  carry  out  the  evaluation  of  the  integral,  nor  does  he  pre- 
sent any  calculations  of  excitation  efficiency.   His  experimental 
measurements  were  concerned  with  the  phase  velocity  and  radial  distri- 
bution of  the  surface  wave  mode. 


2.   MATHEMATICAL  FORMULATION 

The  problem  to  be  investigated  is  illustrated  in  Fig.  1.   An 
infinitely  long,  dielectric  rod  of  radius  b  is  located  such  that  its 
longitudinal  axis  corresponds  to  the  z  axis  of  circular  cylinder 
coordinates  (p,  <f>,    z).   The  rod  is  considered  lossless,  (J  =  0,  with 
a  magnetic  permeability  n      and  permittivity  £,  =  €  €  ,  where  €   is  the 
relative  dielectric  constant.   The  medium  surrounding  the  rod  and 
extending  to  infinity  is  free  space,  with  constants  n   and  €  .   The 
electromagnetic  field  source  is  a  filamentary  ring  of  magnetic  current 
located  at  the  plane  z  =  0.   The  ring  is  of  radius  a,  where  0  <  a  <  b, 
and  is  infinitesimally  small  in  cross  section.   The  source  distribution  is 
represented  as  a  product  of  Dirac  delta  functions  in  the  p  and  z 
coordinates  as  follows: 

K  =^5(P-a)  o(z)  (2.D 

where  (j)   is  a  unit  vector  in  the  (j)  direction,   The  source  distribution 
~K   is  independent  of  (j)   (uniform),  and  is  a  unit  source  such  that 


ff* 


a+A   o+  /> 
0da  =  J  j  5(p-a)  §(z)  •  dz  dp  =  1 

a-A   o-A 


a+A  o+  A 

£(p-a)  dp  =  1,  and    j    £(z)  dz  =  1. 
a-A  o- A 

K  has  the  dimensions  of  volts  per  square  meter.   A  harmonic  time 

dependence  of  e     is  selected  for  the  source. 
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r 

02^ 


The  electromagnetic  field  is  a  solution  of  Maxwell  s  equations. 

-iG)t 
Written  in  differential  form  for  e     time  dependence,  we  have 


VX  E  =  +  i'OfjiH  -  K 

(2.2) 
VX  H  =  -  iCi)€  E. 

Taking  the  curl  of  the  second  equation  and  then  substituting  the  first 
relation  for  VX  E  yields 

_    2    —         — 

-  V  X  V  X  H  +0)  !i6  H  =  -  iO€  K.  (2.3) 

The  only  non-zero  component  of  K  is  the  coefficient  of  the  unit  vector 
<f>.      We  may  write  the  0  component  of  the  vector  eq„  (2.3)  and  obtain 
the  non-homogenous  scalar  equation 


[-  v  x  vx  h]  .  +  coVe ■  H  .  =  -  iO>€  £(p-a)  $(z), 


where  the  (j)  component  of  the  bracketed  term  is  indicated. 

The  magnetic  current  filament  generates  a  field  having  components 
H.,  E  ,  and  E  ,  while  the  components  E  .,  H  ,  and  H  are  equal  to  zero. 
Due  to  the  symmetry  of  the  source,  the  field  is  independent  of  <p   and 
we  note  that  all  partial  derivatives  with  respect  to  (f)  must  be  zero.   It 
is  evident  that  the  source  produces  a  field  which  is  circularly 
symmetric  and  transverse  magnetic  with  respect  to  the  z  axis.   Expanding 
the  bracketed  term  in  cylindrical  coordinates,  we  obtain  the  partial 
differential  equation  relating  H.(p,z)  to  the  source 

3  H,     3H.  32H, 

_2  +  i--£+  (k^  -^H.  +  —=2=  -  loe$(p-a)        S(z)         (2.4) 

dp2        p   2p  p2    *     3z2 


2  2 

where  k      =4)     [it, 

From  E   =  -  -r—z  ~7     XH,    it   follows   that 
i0)€ 


3H 

E      (P,z)    =-r^     o-2 
p     ^  iq€     3z 


(2.5) 


E„    (P.z)    ="^Z    (-^+^H.) 


In  order  to  solve  (2.4)  for  H.,  we  shall  apply  the  method  of  integral 

transforms  to  reduce  (2.4)  to  a  non-homogeneous,  ordinary  differential 

8 
equation.   We  define  the  Fourier  transform  of  H.  (P,z)  as 

+  QD 

h(P,C)  =  J        H0  (P'z>  e_1^Z  dz-  (2°6> 


The  inverse  transform  is  given  by 


+oo 


H4> (p,z)  =  h    J    h  (P>C)  e+iCz  d^  (2-7) 


Assuming  that  the  transform  of  H.  (P,z)  exists,  we  multiply  each  term 

-i£z 
of  (2.4)  by  e  b  and   integrate  over  the  infinite  range  with  respect 

to  z  to  obtain 

A  +  1  dh  +  (k2  _  ^2    1  )  h  =  _  lo€   5(p_a)  (2  8) 

dp    p  dp  P 

where  the  source  variation  £(z)  is  no  longer  present  since 

+co 

S(z)  e"1^2  dz  =  1. 
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We  may  consider  (2.8)  as  an  ordinary  differential  equation  in  which  £ 

is  a  parameter  constant.   The  boundary  conditions  on  h  (P,£,)  are 

obtained  by  taking  the  Fourier  transform  of  the  original  boundary 

conditions  on  H.  (P,z)  and  -r-£.   The  problem  is  one  of  solving  (2.8) 
9  OH 

forh(p,£)  subject  to  the  transformed  boundary  conditions,   H.  (P,z) 
is  then  obtained  by  use  of  the  inversion  integral  (2.7) .   The  integral 
expression  for  H.  (P,z)  is  derived  in  the  following  chapter . 
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3.   SOLUTION  OF  THE  BOUNDARY  VALUE  PROBLEM 


Consider  the  form  of  the  non-homogeneous  differential  equation  (2.8), 
One  could  solve  (2o8)  using  the  Hankel  transform;  however,  the  definition 
of  the  delta  function  allows  one  to  solve  (2.8)  in  a  simpler  manner.  The 
delta  function  is  defined  by  the  relations 


a+A 

S 

a-A 


6(pra)  dp  =  1,  and  6(p-a)  =  0  for  p  /  a. 


Consequently,  for  all  p  other  than  p  =  a,  equation  (2.8)  reduces  to  the 
homogeneous  differential  equation 


d  h   1  dh   ,,2   -2  1  x  ,_ 


(3.1) 


The  delta  function  6(p-a)  implies  a  boundary  condition  which  h(p,£) 
must  satisfy  at  p  =  a.   Multiplying  each  term  of  (2.8)  by  dp  and  inte- 
grating over  the  interval  2A  from  p=a-Atop=a+A,  one  obtains 


dh 

dp 


a+A 


-I  a-A 


,  ,  !<*>«,  +  (k» .  t«, 


h  dp 


a+A 


a+A 


i(^)  dp  =  -1U8  J     S(p-a)  dp 


-A  a-A 

(3.2) 
Assuming  that  h(p,£)  is  continuous  for  all  p,  in  the  limit  as  A— »0 
and  p-»a,  Eq„  (3.2)  reduces  to 


=  -io>8 


(3.3) 


p=a 
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We  have  shown  that  a  continuous  h(p,£>)  which  satisfies  the  homogeneous 
equation  (3.1)  and  whose  first  derivative  is  discontinuous  by  -iwg,  at 
p  =  a  is  a  solution  of  (2.8). 

The  remaining  boundary  conditions  on  h(p,£>)  follow  from  the  boundary 
conditions  imposed  onllAp,z)    and  E  (p,z)  by  Maxwell's  equations.   Refer- 
ring to  Fig.  1,  we  denote  the  cross  sectional  area  of  the  rod  as  region 
I  and  the  space  outside  the  rod  as  region  II.   Since  tangential  H  is 
continuous  at  a  magnetic  current  discontinuity,  we  see  that  Hi(p,z)  is 
continuous  at  p  =  a  for  all  z  including  the  filament  position  z  =  0. 
Since  tangential  E  and  H  are  continuous  across  a  dielectric  boundary, 
we  note  that  H.(p,z)  and  E  (p,z)  must  be  continuous  at  p  =  b  for  all  z. 
The  corresponding  conditions  on  h(p,£)  are: 

(i)   Since  UAp,z)    is  continuous  at  p  =  a  for  all  z, 


pt)  =  J°   H0(p,z)  e-^: 


must  also  be  continuous  at  p  =  a,  hence 


»(P^)]a.A  =  »(P.O]a+A  (3.4) 


(ii)  Similarly,  since  H,(p,z)  is  continuous  at  p  =  b  for  all  z, 


h(P.O]b.A  =  h(p,U]b+A  (3.5) 

(iii)   From  (2.5)  we  write 

l    ^Ed>      l 
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In  region  I,  where  0  -  0  -  b,  the  permittivity  £   =  £  g  .   In  region  II 
where  p  >  b,  £  =  g  .   E  (p,z)  continuous  at  p   =   b  for  all  z  requires 


!5*  i 


=  8. 


?5*   i 

ep  +  p  H0 


b-A 


b+A 


Taking  the  transform,  we  obtain 


dh        1    ,. 
dp+ph 


=  e 


b-A 


dh        1   . 

ap+Ph 


b+A 


which  reduces  to 


Sr  [f]    -  [§]    ♦  [|  >]   C8r-X)  -  o 

L^J  b+A    L^J  b-A    LH  J  b+A 


(3.6) 


by  reason  of  (3.5) 


In  order  to  determine  H*(p,z)  in  region  II,  we  must  solve  Eq.  (2.8) 
for  the  corresponding  h(p,£>).   Inside  the  dielectric  rod,  for  all  p 
except  p  =  a,  (2.8)  becomes 


dp 


h   1  dh   ,,2 
2  +  p  dp  +  (kl 


^   -   i)  h  -  0 


(3.7) 


where  k,  z  =  CO' 


^O&i 


Outside  the  dielectric  rod  and  for  all  p,  (2.8)  becomes 


dh   ldh    .2    ,2    1  N  . 

~2  +  p  dp  +  (k0   -  &   "  "2>  h  =  ° 
dp    ^  ^  p 


(3.8) 


where  kQ2  =0)  fi0&(). 
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We  proceed  by  writing  general  solutions  h(p,£>)  of  (3.7)  for  the 

regions  0  -  p  -  a,  and  a  -  p  -   b.   The  discontinuity  at  p  =  a  appears 
in  the  boundary  condition  (3.3).   A  general  solution  of  (3.8)  is  written 
which  yields  h(p,£)  for  the  region  p  -  b.   The  three  solutions  possess 
six  arbitrary  constants.   The  constants  are  determined  from  the  require- 
ments that  the  field  be  finite  at  p  =0,  that  the  field  be  regular  at 
infinity  and  that  it  satisfy  the  radiation  condition  and  the  four  boundary 
conditions  (3.3),  (3.4),  (3.5),  and  (3.6).   Choosing  the  solution  h(p,£) 
for  region  II,  the  inversion  integral  (2.7)  must  then  be  evaluated  to 
yield  H.(p,z)  for  p  >  b. 

Equation  (3.7)  is  recognized  as  a  form  of  Bessel's  differential 
equation.   It  has  the  general  solution 


h(p,£)  =  A  J^T^p)  +  P  y^-tJ^p) 

where  A  and  P  are  arbitrary  constants,  J  (li  p)  and  Y  C\)   p)  are  Bessel 
functions  of  the  first  and  second  kind,  respectively,  and 

Consider  the  case  0  -  p  -  a.   Since  Hj(p,z)  and,  therefore,  h(p,£) 
must  be  finite  at  p  =  Q,  we  determine  that  P  =  0,  since  Y  (V  p)  is 
unbounded  as  p— >0. 
Hence 

h(p,£)  =  A  Jx  Ci^p)  (3,9) 

where 


0  6  p  £ 


^=[kl-^] 


2_l/2 


M 


For  the  region  a  -  p  -  b,  the  Y   function  must  be  included,  and  we  have 


where  a  -  p   -  b; 


h(p,t,)  =  B  J1Cl)1P)  +  C  Y±    {])  p)  (3.10) 


P  1/2 


For  the  region  b  -  p  we  choose  to  write  the  general  solution  in  terms  of 
Hankel  functions  of  the  first  and  second  kind 

h(p,;>  =  d  h1(1)  ay))  +  q  h1(2)  g)^) 

where 

^o  -  K2  -  tf* 

It  will  be  seen  later  that,  in  general , "U   is  complex;  therefore  we 
define  the  argument  ofi)   to  be  0  —  arg  l)     -  IT.     As  P-=»oo,  H     ^cp^ 

(2) 

vanishes  and  H     ^c\P}   is  unbounded  for  the  defined  argument  ofVft; 


therefore  Q  =  0  and 


h(p,t,)  =  D  H1U)  (VoP)  (3.11) 


where 

b  ^  p 


^o  "  tko  "  tl 


0  -  arg-u>Q  -  7L 


In  order  to  solve  for  the  arbitrary  constants  A,  B,  C,  and  D,  we 
apply  the  boundary  conditions  (3.3),  (3.4),  (3.5),  and  (3.6)  to  the 
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appropriate  solutions  (3.9),  (3.10),  and  (3.11).   Omitting  the  details, 
the  following  four  equations  are  obtained: 


-  (vxb)  A  Jx'  (vxa)  +  (v^)  B  J^  (v^)  +  (v^)  C  Y^  (v^)  =-x^±b 


A  J1(via)  -  B  J^i^a)  -  C  Y  (v^a)  =  0 


B  J^v^b)  +  C  Y1(vib)  -  D  H^V^b)  =  0 


>ib)  B  j'  (vib)  +  (vib)  C  Y^  (v^)  -  D[£r(vQb)  H1/(vQb)  + 

+(6r-D  Vvob)]  =  0 

(3.12) 

This  system  of  equations  was  written  in  matrix  form  and  solved  by 
finding  the  inverse  of  the  coefficient  matrix.   Details  are  given  in 
Appendix  A.   The  constant  D  was  determined  to  be 


Jx(va) 

D  =  1UV (I) (I) 

(vxb)  J0(vib)  H/  ;(vQb)  -£r(vob)  J^b)  H0C  ;(vQb) 

(3.13) 

Substituting  (3.13)  into  (3.11)  we  obtain  h(p,?,)  for  the  region  p-  b. 

The  solution  h(p,£>)  may  then  be  substituted  into  the  inverse  trans-form 

(2.7)  to  yield  the  integral  expression  for  Hj(P,z)  in  region  II.   This 

result   is  .    . 

icj^a         p+0°  J^a)   H^V^)    e+i^z  d£ 

VP'Z>  ^      (v^J^b)   Hl(1)(v0b)   -Sr(v0b)   J^b)    Hoa)(v0b) 

(3.14) 
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where  p  ^  b  , 

vi  =  Lki   - 1  J     j  V0  -  Lk0   - 1  J   • 

The  evaluation  of  (3.14)  is  the  subject  of  Chapter  5. 
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4.   SINGULARITIES  OF  THE  INTEGRANT*—  SOLUTION  OF  THE  MODE  EQUATION 

The  integral  for  H .  (p,  z)  is  treated  as  a  contour  integral  in 
the  complex  t,   plane  where  the  path  of  integration  is  along  the  real 
axis  from  -co  to  +co  .   In  order  to  carry  out  the  integration,  it  is 
first  necessary  to  determine  the  singularities  of  the  integrand.   A 
detailed  analysis  of  the  singularities  is  presented  in  this  chapter. 
The  expression  for  H.  (P,z)  is 

+i0)€ia    +^    J^a)  Hl(1)  (VQP)  e+iCz  dC 


H.(p,z)  =  —  -        ,    t^t ^p- 

211  ^  (Vlb>      Jo<Vlb>     *!  <Vob>      -     MV„b>VVlb>0  (V     b> 

-oo     1     ol    1      o      rollo     o 

(4.1) 

4.1  Branch  Points 

Consider  the  variable  v,  which  appears  in  the  arguments  of  J  and 
1 o 

I   2     2 
J  in  (4.1).   Since  v  =  k   -  £  ,  then  v  is  multiple-valued  in  any 

neighborhood  of  I,   =  +  k  .   However,  if  one  considers  the  power  series 

expansions  of  J  and  J  ,  one  sees  that  the  integrand  is  an  even  function 

of  v  so  that  the  points  t,   =  +:  k   are  not  actually  branch  points.   The 

I  2  2 

Hankel  function  arguments  contain  the  variable  v  =  k   -  t,    .   The 

Hankel  function  has  a  logarithmic  singularity  at  v  =0,  and  so  is 

multiple-valued  in  any  neighborhood  of  t,  -   +  k  .   The  points  £  =  +  k 

are,  therefqre,  branch  points  of  the  integrand.   We  select  branch 

cuts  in  the  t,  plane  as  shown  in  Fig.  2.   The  points  +  k  and  +  k  lie 

on  the  real  axis  with  I  k,  I  >  '  k  I  since  € ,  >  €  .   The  branch  cuts 

I  1|  '  I  o  I        1     o 

are  defined  as  t  =  +   k  f  i  In  1. 
—  o 
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4  2   Poles 

Before  we  investigate  the  integrand  of  (4  1)  for  poles,  it  will 
be  advantageous  to  consider  the  significance  of  a  pole  in  the  physical 
problem.   Assume  that  the  point  t,   =  £   is  a  pole  of  the  integrand. 
When  (4.1)  is  evaluated  as  a  contour  integral,  a  residue  contribution 
at  £  =  4  must  be  included.   Thus,  Hj(P,z)  will  have  a  term  of  the  form 

nx        it  z 

Ah/A;  (v  p)  e   °  (4.2) 

1      o 


where  A   =  27Ti   Residue 


*-V 


2    2 
Since  we  assumed  a  lossless  dielectric  rod  by  taking  k  =  03  \jl€   real,  we 

will  expect  (4.2)  to  represent  a  wave  propagating  without  attenuation, 

which  implies  that  t   is  real.   If  t   is  real  and  It  I  >  |k  |  , 

o  o  I  o  |     I  o  I 

then  v   is  pure  imaginary  and  H     (v  P)  reduces  to  the  K.  Bessel 

function  which  decays  exponentially  with  increasing  argument.   Under 

these  conditions,  (4.2)  represents  a  surface  wave  of  amplitude  A 

having  a  radial  dependence  of  H     (v  P)  which  propagates  along  the 

-i(CJt-£oz)    l  ° 

rod  according  to  e  .We  shall  now  determine  if  there  are 

values  of  t,   real  which  cause  the  denominator  of  (4.1)  to  vanish. 

Equating  the  denominator  of  (4.1)  to  zero  yields 

J  (v.b)  H  (1)(v  b) 

1     o 

For  convenience  of  notation,  let  v   b  =  X,  and  v  b  =  X  ;  then  (4.3)  may 

11  o  o 


be  written 


J    (X,)  H    (1)(X   ) 

v        ol      =     g  o  o 

1  W  r      °Hl(1)(x) 

1  o 


Recalling  that 


vl = 

k2  - 1 

o 

K2  -  <? 

2  2 

where  k   >  k   >  0, 
1     o  y 

(i)   Consider  the  case  of  £  real  and  It,  |  >  Ik  I  >  Ik  I 


Then 


v,  and  v  ,  and  therefore  X,  and  X  will  be  pure  imaginary „   Denote 
1      o  1      o 

X,  =  ip  and  X  =  i£ ,  where  p  and  £  are  real;  then  (4.4)  becomes 
1  o 

J  (ip)  H  (1)(i£) 


(iP)  tW-  <L  ("i£)  — ; 


which  may  be  reduced  to 


Vip)     r       Hl(1)(i4) 


i  (p)  k  (i) 


where  I   and  K  are  modified  Bessel  functions  of  the  first  and  second 
n      n 

kinds,  respectively.   The  left  hand  side  of  (4,5)  is  always  negative 
for  real  p  while  the  right  side  is  always  positive  for  real  |,  so 
there  are  no  values  (p,G)  for  which  (4,5)  can  be  satisifed.   We 
conclude  that  no  poles  can  occur  on  the  real  axis  for  I  £  I  y   Ik  I. 
(ii)   Consider  t,   real  and  £|  <  |k  I  /  Ik  1  ;  then  v 


and  v  . 
o 


and  therefore  X,  and  X  are  pure  real.   From  Eq.  (4.4)  we  write 
1      o 

J   (X>  H  (1)  (X  )  H   (X  ) 

1      o  00 

where  H*  denotes  the  derivative  of  H 
o  o 
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Using  the  definition  of  the  Hankel  functions,  (4.6)  may  be  written 

^V  _  .  £   x   W  -  |  W 

Al  Jn  (Xn  )  ~     r   o   J  /  (X  )  +  i  Y  /  (X  )  "  K*'" 

11  o   o       o   o 

Since  X   is  real,  the  left  side  of  (4.7)  is  always  real,  while  the 

right  side  is  obviously  complex  for  real  X  .   Is  it  possible  for  the 

o 

right  side  to  be  real  for  some  X  ?  The  right  side  can  be  real  only 
if  the  argument  of  the  denominator  is  equal  to  that  of  the  numerator 
by  a  multiple  of  nTT. 
Let 


a  =  arg  [j^)  +  i  Y^)] 
P  =  arg   [Jo(Xo)  +  i  Yo'(Xo)] 


Assume  |3  =  a  +  nlT, 

where  n  =  0  +  1 ,    +  2 , 


then 


since 


a        4.         ,  «\         tan  a       +   tan  nff 

tan  p   =  tan    (  a+nTT)    =  - — — — — =-=  tan  a 

r  1   -   tan  a  tan  nTT 


Y      (X    )  Yj{Xn) 

tan  a  =  and   tan  [3 


We   obtain 


J      (X    )  J  '   (X   ) 

o        o'  o        o 


Y    (X    )         Y  >  (X   ) 
o     o  o        o 


J    (X   )         J  /  (X   ) 
o      o  o        o 


J    (X   )    Y  '     (X   )    -   Y    (X   )    J  I  (X   )    =  0.  (4.8) 

oooooooo 

However,    the  left   side  of    (4.8)    is   the  Wronskian  of  J   ,    Y   ,    which 

o   o 

is  non-vanishing  for  finite  X  : 

o 


W   £J   ,Y  ?    =  J    (X   )   Y'    (X   )    -Y   (X   )J    \x   )    =  -z§-  *  0 
t  o      OJ  oooo  oooO  7TX 
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We  conclude  that  the  right  side  of  (4„7)  must  always  be  complex. 

Since  the  left  side  of  (4.7)  can  vanish,  is  it  possible  for  the  right 

side  to  be  zero?  For  the  right  side  to  be  zero,  J  (X  )  and  Y  (X  ) 

o   o       o   o 

9 
must  be  zero  simultaneously,   Watson  proves  that  the  positive  zeros  of 

two  distinct  cylinder  functions  of  the  same  order  are  interlaced; 

therefore,  J  (X  )  and  Y  (X  )  are  never  zero  simultaneously.   The 

same  holds  for  -X  since  J  (-X  )  =  J  (X  )  and  Y  (-X  )  =  Y  (X  ) . 
o        o   o     o  o       o   o     o  o 

Since  the  right  side  of  (4.7)  is  always  complex  and  never  vanishes, 

there  are  no  real  values  X.  and  X  for  which  (4,6)  can  be  satisfied. 

1      o 


We  conclude  no  poles  exist  for  t,   real  and  I  £  <  Ik 


(iii)  Consider  Z,   real  and  Ik  I  <  |£|  /  Ik,  I  ;  then  y  is  real,v 

is  pure  imaginary,  and  therefore  X,  is  real  and  X   is  pure  imaginary. 

1  o 

Let  X  =  ib ,  where  b  is  positive  and  real,   Equation  (4.4)  becomes 

i(V 


J  (X,)  H  (1)(i£)        K  (£) 

=     €  (-i£)  -°-- =  6  e  ^TT  .  (4.9) 


jo(xx)       ko(£) 

A  graph  of  the  functions  -X     .   .  and  €  £    g   is  presented  in  Fig.  3. 

1  Jn  (X  )  r   Kn  (5 ) 


Jo  (X1} 

The  function  -  X     .      >   is  discontinuous  each  time  J, (X, )  vanishes 

1   1 


lv  1 


and  so  an  infinity  of  branches  occur  for  ever  increasing  X  .   It 

passes  through  zero  each  time  J  (X,  )  has  a  zero.   The  function 

o  1 

KQ(b) 

€  £  v   (t\    approaches  zero  for  b — »0  and  increases  positively  for 

increasing  £.   It  is  evident  in  Fig.  3  that  for  every  finite  b,  an 
infinity  of  values  X  exist  which  satisfy  Eq«  (4.9).   In  order  to 
obtain  a  finite  number  of  unique  solutions  of  the  mode  equation, 
a  second  relation  between  X1  and  £  is  needed. 


2;j 
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It   follows   from   the  definitions 


A2"2 


,     2         r2 

v     = .  k       -  t, 

o       \o 


where  £  is  identically  the  same  in  v     and  v  . 


,2 

Elimination  of  ?,   yields 


2 
2  2  2  2  27T 

v/  -  V /   =  k_-       -  k         =    (J-)       (€      -    1) 
1  o  1  o  A.  r 

o 


2 

Multiplying  by   b      yields 


(vlb)2  -  (vob)2  =  C2^)2  (€r-l) 
o 


o 


Si 


nee  X  =  i£ ,  we  obtain  the  second  relation  between  X  and  %  ,    which  is 

oqq 
Xr   +  i      =  R  (4.10) 


where 


,27Tb 


)JVi   =kobjvr 


*  HP)   € 


The  graphical  solution  of  the  mode  Eq.  (4.9)  is  illustrated  in  Fig.  4, 
which  includes  curves  of  X  as  a  function  of  §  obtained  from  Fig.  3, 
and  the  relation  (4.10).   Note  that  (4.10)  is  the  equation  of  a  circle 
of  radius  R  with  center  at  the  origin. 

The  multiple-valued  solutions  X  for  every  £  as  discussed  with 
Fig.  3  are  evident  in  Fig.  4.   The  first  branch  of  X  starts  at  X   =  2.405 
for  t    =   0  and  approaches  X  =3.83  asymptotically  with  increasing  |. 
Xx    =   2.405  is  the  first  zero  of  J^)  and  X  =  3,83  is  the  first  zero 
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FIGURE  4.      GRAPH  OF  X     =    f(^)    and   X  2  +   £2  =   R 
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of  J  (X  ) ,   The  second  branch  commences  at  5.52  and  approaches  7,02, 


which  corresponds  to  the  second  zeros  of  J  (X  )  and  J  (X  ) ,  respectively 
An  infinity  of  branches  is  thus  established 


Consider  the  case  R  <  2.405.   Since  X±    =  \|  R  -  £2  ^   R,  Fig,  3 
shows  that  Eq.  (4.9)  cannot  be  satisfied  for  X  <  2.405,  that  is,, 
no  solution  of  the  mode  equation  exists  for  R  <  2.405.   This  is  also 
evident  in  Fig.  4,   Physically;  this  means  that  the  dielectric  rod 
waveguide  is  below  cutoff  and  cannot  propagate  a  surface  wave  of 
the  transverse  magnetic  type. 

Referring  to  Fig.  4,  if  2.405  <  R  <5.52,  a  single,  unique  solution 

(£,    X  )  of  (4.9)  and  (4.10)  results.   Thus,  (4.1)  has  poles  at  t,   =  +  t,    , 

where  t,     may  be  determined  from  the  solution  (£  ,  X_  )  „   The  poles  occur 

on  the  real  t,   axis  in  the  region  |  k  J  /    \t,\  •  Ik  I  .   We  see  that  when 

R  =  Pr— ) . |€   -1  is  restricted  to  the  range  2.405  <  R  <  5.52,  the 

dielectric  rod  propagates  a  single,  surface  wave  which  is  the  lowest 

order,  circularly  symmetric,  transverse  magnetic  mode.   It  is  known 

i£oZ  -i£oz 

as  the  E   mode  and  propagates  as  e     for  z  -7  o  and  e      for  z  <  o. 

As  R  increases, the  number  of  solutions  (£ ,  X  )  of  (4.9)  and  (4.10) 

also  increases,  and  the  dielectric  rod  supports  an  ever  increasing 

number  of  higher  order  TM  modes.   The  number  of  poles  occurring 

between  k»  and  k  increases  with  R,  as  is  shown  in  Fig.  5.   It  is 

interesting  to  note  the  behavior  of  k  ,  k  ,  and  the  poles  as  R  increases. 

From  the  definition  of  v  we  obtain 
o 


.  k  2  +  (|)2  (4.11) 

\|  o     b 
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Consider  R  just  greater  than  2.405;  then  X  ft*  2.405  and  b  is  a  very 
small  positive  quantity.   Equation  (4.11)  shows  £   is  just  greater 
than  k    ,  that  is,  the  poles  +  £  first  appear  in  the  neighborhood 
of  the  branch  points  +  k  with  \t,\     ~7       k   ,   Referring  to  Fig.  4, 
as  R  increases,  the  s  of  any  solution  increases,  so  that   £,  corresponding 
to  (£ ,  X  )  also  increases.   Thus,  the  poles  +  £  move  away  from  the 
origin.   Each  new  pair  of  poles  +  t,      first  appear  near  t,   =  +   k  ,   Now 


27T         1 I 

k   =  (r— ) ,  k,  =  J6   k  ,  and  R  =  k^b  €-  -  1,  so  we  may  think  of  k   and  k 
o    \  1   M  r  0  6  \J  r  '  1       0 

as  proportional  to  R.   The  distance  from  k   to  k  ,  (k  -k  )  =  k   (  J€~-1) 

is  also  proportional  to  R.   As  R  increases,  k  and  k  move  away  from 

the  origin,  the  interval  (k  -k  )  increases,  and  the  number  of  poles 

between  k  and  k  increases. 

4.3   Solution  of  the  Mode  Equation 

The  object  of  this  work  is  to  determine  the  efficiency  of  a  magnetic 

current  filament  in  exciting  the  E  _  mode  on  a  dielectric  rod.   For 

single  mode  operation,  R  must  be  restricted  to  the  range  2.405  <  R  <  5.52, 

Since  R  =  k  b  I  €  -1 ,  this  restriction  on  R  is  actually  a  limitation  on 

k  b  and  €  ,  the  physical  constants  which  describe  the  dielectric  rod 
0       r 

as  a  waveguide.   Specific  solutions  of  (4.9)  and  (4.10)  are  required 
in  order  to  calculate  the  power  transmitted  in  the  surface  wave  and, 
ultimately,  efficiency.   Equations  (4.9)  and  (4.10)  were  solved  for 
the  case  €   =  2.56,  a  relative  dielectric  constant  typical  for 
polystyrene.   Six  solutions  (5 ,  X  )  were  obtained  for  selected  values 
of  k  b  equal  to  2.2,  2.6,  3.0,  3.4,  3.8,  and  4.2.   Graphical  solution 
of  (4.9)  and  (4.10)  as  illustrated  in  Fig.  4  was  considered  too 
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inaccurate  for  the  purposes  of  numerical  computation  and  so  (4.9) 

and  (4.10)  were  solved,  utilizing  the  University  of  Illinois  digital 

computer.   The  solutions  (£ ,    X±)   are  accurate  to  four  decimal  places. 

Details  of  the  method  used  to  solve  the  equations  on  the  computer 

are  given  in  Appendix  B.   Table  I  presents  i,    X  ,  and  the  ratio 

W\   corresponding  to  each  value  of  k  b.   The  ratio  X  |X  follows 
s '  0  g|  o 

from  the  definition  of  v  . 


X  |X_  = 
g 


"  M^©3 


(4,12) 


TABLE  I 


V 


2.2 
2.6 
3.0 
3.4 
3.8 
4.2 


5 

Xl 

x  rxn 
g|  o 

5603 

2.6901 

.9691 

1.2329 

3.0043 

.9036 

1.9353 

3.2086 

.8403 

2.6269 

3.3366 

.7913 

3.2905 

3 . 4204 

.7560 

3 . 9264 

3.4788 

.7305 
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5,   EVALUATION  OF  THE  CONTOUR  INTEGRAL 

For  the  purposes  of  the  following  discussion,  we  repeat  the  integral 
solution  for  H-(p,z)  outside  the  dielectric  rod,  which  is 

+iG)&a    p              J.(va)  H(1)  (vnp)e+i^>z 
VP'Z)  "  IT  J  -IT) " ' (1)     di; 

c0(vib)  Jo(V»  Hi      <vob)  "  W>  Ji(V>  Ho  V> 

(5.1) 

where  C   is  the  contour  along  the  real  axis  with  appropriate  indentations 
o 

+  +  * 

at  the  branch  points  -  k  and  poles  -  b„,  as  shown  in  Fig.  2.   Expression 
o  0 

(5.1)  is  evaluated  by  transforming  from  the  £  plane  to  the  complex  -r  plane 

and  employing  the  saddle  point  method  of  integration  to  obtain  an  asymp- 

10,11 
totic  value  of  the  integral.      The  saddle  point  method  is  applicable 

only  when  certain  parameters  in  the  integrand  are  relatively  large.  It  is 
used  to  obtain  an  approximation  to  the  far  zone  radiation  field.   The  sur- 
face wave  field  follows  from  the  residue  contribution  at  the  pole.   In  the 
T  plane  the  contour  of  integration  is  deformed  until  it  passes  through  the 
saddle  point  of  the  integrand.   The  integrand  contains  an  exponential 
function,  and  when  the  deformed  contour  corresponds  to  the  path  of  wsteep- 
est  descent   the  exponential  function  decreases  very  rapidly  away  from  the 
saddle  point.   The  integrand  vanishes  at  the  ends  of  the  path,  and  so  the 
principal  contribution  to  the  integral  comes  from  evaluation  along  the  path 
of  steepest  descent  in  the  very  near  vicinity  of  the  saddle  point. 
5.1   The  Radiation  Field 

In  order  to  apply  the  saddle  point  method  of  integration,  we  introduce 
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the  transformation  of  variable 

£>  =  k   sin  t 
o 

where 

t  =  4«  +  iT|  (5.2) 

Now  t  =  sin   (£|k  )  is  a  multiple-valued  function  of  £,  and  the  region  of 
integration  in  the  C  plane  transforms  into  a  strip  in  the  T  plane  bound 
by  two  curved  lines  which  correspond  to  the  branch  cuts  in  the  £,  plane. 
They  are  defined  by  the  equations 

C  sin  ^  coshTj  =1  C  sin  ^   cosh  T\  =  ~1 

1  L  ^  I  L 

(^  cos  +  sinh  7]   =  ^-2  >  0  L  cos  ^   sinh  r\  =  ™-^  ^    0  (  5 .  3) 


o  o 


The  contour  of  integration  C  in  the  ^  plane,  which  is  defined  as  Im  %,   =  0 

with  necessary  detours,  transforms  to  the  path  C   in  the  t  plane,  as  shown 

in  Fig.  6.   It  follows  from  the  definition  Im  £  =  0  which  corresponds  to 

k  cos  4*  sinh  71  =  0  in  the  t  plane  while  maintaining  C  within  the  curved 
o  '  1 

strip.   The  direction  of  C  follows  from  consideration  of  the  real  part 

of  £  along  C_,  where  Re  £  =  k   sin  41  coshTh   The  branch  points  £  =  -  k 
U  o  v  o 

transform  to  the  points  t  =  -  7T/2  while  the  images  of  the  poles  -  fe  are 

the  points  T   =  +7T/2  -  i  cosh   C  /k  and  -f  =  -7T  2   +  i  cosh   £»  /k  . 
o     '  o  /  o       o     '  o'  o 

Substituting  (5.2)  into  the  expression  for  v  yields  v  =  Jk2  -  C2  - 

-  k   cos  1"  „   The  choice  of  sign  here  is  determined  by  the  previous 

definition  of  arg  v   in  the  £>  plane.   On  the  contour  C_  ,  %,   is  real  and 
o  U 

when  |£>|  >  k  I  we  obtain       v      =    j  v    I  g  ±  *•  V2 .   Since  we  chose  to 

define  0  <  arg  v  <  77",  the  negative  sign  is  excluded  and  we  have 

V0  =  Vn  e+  1     ■   In  orde^  that  corresponding  values  of  f  on  C  yield 
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the  same  argument  for  v    the  positive  sign  is  required;,  thus  we  establish 


vn   =  +  k  cos  t.  (5.4) 

0       o 


2    Y2        +  2 

The  function  v,  becomes  v      =   k   -  fe   =  -  k   ,  £   -  sin  t.   Since  the 
1  1    \|  1  0  \j^r 

integrand  of  (5.1)  is  an  even  function  of  v    ,  either  sign  may  be  used 
and  we  choose 


2 
-  sin  t  =  k  w(t)  (5.5) 

o  \J  r  o 

where  for  convenience  of  notation  we  let  w,  a  function  of  t,  represent 

the  radical,  w  =   £  -  sin2  T.   Substitution  of  (5.2),  (5.4),  and  (5,5) 

into  (5.1)  with  d£  =  k  cos  t  dT  yields  the  transformed  integral  in  the  t 
o 

plane  which  is 


ico£  a  n  J,  (k  a  w)  H  (1)(k  p   cos  t)  elk°Z  sinT  cos  T  dr 

rr/^\  1  lO  1  O 

H  .  (P  ,  Z)   =  T-T pj 

v  27Tb       U     w  J  (k  b  w)   H,         (k  b  cos  t)-  £      cos  t  J   (k  bw)  Bt>    (kb  cost) 

Do  lo  r  louo 

Cl 

(5.6) 


If  p  is  large  and  k  p  cos  T  =f   0,  the  Hankel  function  may  be  replaced  with 
its  asymptotic  representation,   Then  the  Hankel  function  and  exponential 
in  (5.6)  combine  in  the  form 

Hl(1)(k  p  cos  T)eik°Z  Sin  T~  e"1?^ 2 >     eiko(P  cos  T  +  z  sin  T> 

or  \  7Tk  p  cos  t  / 

(5.7) 
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Consider  the  spherical  coordinate  system  illustrated  in  Fig.  7,  where  the 


FIGURE  7   CIRCULAR  CYLINDRICAL  COORDINATE  SYSTEM  (p,0,z)  AND  SPHERICAL 
COORDINATE  SYSTEM  (r,0,9) 


polar  angle  Q  is  measured  from  the  plane  z  =  0.   In  this  coordinate  system 
we  note  that  p  =  r  cos  9  and  z  =  r  sin  9.   Substituting  for  p  and  z  in  the 
right  side  of  (5.7)  we  obtain 


Hn(1)(k  p  cosT) 
1     o 


37T 


ik0z  sinT  ^  ft-i-£-  p        2   ; 

r  cos  6  cos  f- 


Lit  k 


V2 


ikQr  cos(-f-0) 


(5.8) 
We  substitute  (5.8)  into  (5.6)  to  obtain  the  following  approximation  for 
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for  H.  (r,9)  which  holds  for  r  large  and  k  r  cos  9  #  0  „ 

i/i 

iooe  a  37T       I  „  ,      ^ 

H*<r,0)    =  — J-  e"1^  F(r)  [      *   C°S   TQ  ]       e1   kor   ^(t-9)    dT 

0      '  27Tb  I  L  77k   r   cos   9  J 


where 

J    (k     a  w) 

F(t)    =  — 2 

w  j_ (k  b  w)   H  v    '    (k  b  cos  t)    -  g      cos  t  jr   (k  b  w)   H_v      (k  b  cos  t) 
Oo  1  o  r  lo  Oo 


Integral  (5.9)  is  readily  evaluated  by  the  saddle  point  method.   The  saddle 
point  is  found  by  setting 

4~   cos  (t  -  9)  =0  (5.10) 

dT 

which  gives 

t  =  e.  (5.11) 

Note  that  the  saddle  point  corresponds  to  the  polar  angle  9.   The  path  of 
steepest  descent  is  defined  by 

Im  [i  cos  (t  -  e)J  =  cos  (41  -  9)  coshT]  =  1, 


which  may  be  written 


sin  (i|j  +  TT/2  -   9)  coshTj  =  1.  (5.12) 


Thus,  the  steepest  descent  path  C  has  the  same  form  as  the  curve  bounding 

s 

the  right  side  of  the  strip,  but  it  is  displaced  to  the  left  by  the  amount 


(7772  -9),  as  shown  in  Fig.  6.   One  may  verify  that  on  C 


dTl 
d$ 


=  -1  at  T  =  e. 
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so  that  the  steepest  descent  path  crosses  the  real  axis  at  the  angle  -7T/4. 
We  deform  contour  C^  into  the  path  C  which  starts  at  0  -  7T/2  +  i  oo  and 
ends  at  e  +  7T/2  -  i  oo  .   In  order  for  (5.9)  to  converge  as  T|  approaches  -  oo 
on  C  ,  it  is  necessary  that 

Re  [i  cos  (T  -0)3  =  sin  (^  -  0)  sinhT)  <  0.         (5.13) 
This  condition  is  satisfied  in  the  region 

T)>0   for  -7T+0<4J<© 
and 

7]<0    for   0    <  4<    <  IT  J  (5.14) 

The  shaded  area  of  Fig.  6  is  the  region  of  convergence  and  we  see  that  the 

integrand  of  (5.9)  does  vanish  on  C  as  T)  approaches  -  oo  . 

The  relation  of  the  path  of  steepest  descent  and  the  polern  will 

now  be  considered.   Evaluation  of  (5.9)  along  C  yields  the  total  field 

s 

in  the  region  outside  the  dielectric  rod.   Since  the  saddle  point  is 

equal  to  the  polar  angle  0,  as  0  increases  towards  7T/2,  the  steepest 

descent  path  shifts  to  the  right  in  Fig.  6.   Ultimately,  a  residue  term 

must  be  included  because  C  will  cross  to  the  right  of  the  pole  T«.   There- 

s  0 

fore,  the  total  field  outside  the  rod,  which  is  obtained  as  0  takes  on 
values  0  ^0  <  tt/2s   may  be  considered  as  the  sum  of  the  residue  at  the  pole 
and  the  steepest  descent  contribution  at  the  saddle  point,  that  is,  the 
surface  wave  and  the  radiation  field . 

Due  to  the  rapid  decay  of  the  exponential  term  away  from  the  saddle 

point,  it  is  only  necessary  to  integrate  over  a  short  segment  of  C  in  the 

s 

near  vicinity  of  Y=  6. 
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We  proceed  by   setting 


-T   -   9  =  a  e"1* 
.7T 
dX  =  do-  e~17^ 


(5.15) 


where  a  is  the  distance  along  C  measured  from  the  saddle  point.   Consider- 

s 

ing  \  f  -   -9|  small,  we  may  approximate 

2         2     7T 
cos  (T-  ©)  ,*  1  -  {rf~2B)-   =  1  -  "Y  e_i2  (5.16) 


Except  for  the  exponential  term,  the  integrand  of  (5.9)  is  slowly  varying 
about  T1  =  @  and  may  be  assumed  as  constant.   We  substitute  (5.15)  and  (5.16) 
into  (5.9)  to  obtain 


Vr'S)  =  «b 


L.  [_|_]   F(e)  elko 
L7T  k  rJ 


+A 


ik^r 


k  ro^ 
o 


dcr, 


(5.17) 


Since  r  is  very  large,  the  following  approximation  is  valid 


J    e-V-d(r«     /     e-V-do.  =  [|i] 


*6 


(5.18) 


Substituting  (5.18)  into  (5.17)  and  noting  that  W£   =  k  £ 

1    or 

the  far  zone  radiation  field,  which  is 


^o 


,  we  obtain 


yr,6)  = 


.  OJ    7Tb        r 
\^o 


(5.19) 


where 


F(0) 


J  (k  a  w) 
1   o 


wl(k  bw)R       (k  b  cos  0)-£   cos  0  J,  (k  bw)    H  (k  b  cos  ©) 

Oolo  r  loOo 


IH 


and 


S   -  sin  9 
r 


It  follows  from  Maxwell  s  equations  that 


Ee  (r'e)  =  -"-S->  p(9)  I 


(5.20) 


5.2  The  Surface  Wave  Field 

The  surface  wave  propagating  in  the  positive  z  direction  is  obtained 

by  evaluating  the  residue  of  (5.1)  at  the  pole  ?>  .   Denote  the  denominator 

of  (5.1)  by  G(£)  so  that  (5.1)  may  be  written 


yp.«>  -w 


f^  (1)         iCz 

ico^a  J^a)  H^  '    (vQp)  elb 


27T 


G  (&) 


d£ 


(5.21) 


where 


G(t)  =  (v1b)  J0(vali)  H1(1)(v()b)  -£r(vQb)  J1(v1b)  H0(1)(v0b) 


The  integral  (5.21)  is  the  solution  for  the  total  field  H,(p,z)  in  Region  II 
where  p  >  b.   Let  H   (Pi>z)  represent  the  surface  wave  portion  of  the  field 
in  Region  II.   It  is  given  by 


H/X   (P>z)  =  27ri  Residue 


(1) 


j^a  Jj^a)  Hxv    (vQp)  e 
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G(&) 


(5.22) 


t-6. 


where  z  >  0. 

The  pole  £   follows  from  the  solution  (I  ,  X  )  of  the  mode  equation  G(£)  =  0, 


Since  G'(&o)  =^  G(&) 


£=  0,  C   is  a  simple  pole  and  the  residue  of 


&  »&, 


of  the  integrand  at  £   is  given  by 


Res 


iut..   Ji(5VH(1),i^)eV 

1  b  b .  (5.23) 


0'<V 

Using  Bessel    function  recurrence   relations,    G   (£ft)   was   determined  as 

2   ,       2j„  (B/g)2 tx        r\         «,     ,  I 


(V  -iW^VV  Ki<e>  +t(T>  +&r(r^  W  Ko(l)  + 


1  1 


+   (6     -    1)    J,     (X.)    K.     (I)    >  (5.24) 

r  1       1       1  / 

Using  the  relation  H  (i  £  £)    =  -  ~  K   (§^),    it   follows   from   (5.22), 

(5.23),    and   (5.24)    that 


H0II(P,z)    =An    Kx    (£^)    e^oz  (5.25) 


where  p  >    b,    z  >  0 


^  (=>  j,  <f  v 


V{2   (xf  J0(X1)K1(I)    +    l>r(4>+(r)]   J0(Xl'   *><e>*<6,-»   WV6)} 


with 


and 


&0b  =  J  ( V)2  +  42 


R2=Xl2+|2. 


It   is   interesting  to  note  that   the  surface  wave   solution  (5,25) 


exist   since  £>     is  real.      This   situation  results   from  the  simplifying 


40 
assumption  that  the  dielectric  rod  is  lossless.   If  the  finite  loss  of  the 
rod  is  included,  the  constant  k  2  =  d)2  u      £,,(1  +  i^]?")  is  complex  since 
<x  f  0,      It  follows  that  the  poles  -  b0 ,  which  are  defined  by  the  solution 
of  the  mode  equation  for  complex  arguments,  no  longer  lie  on  the  real  axis. 
Instead,  they  appear  in  the  first  and  third  quadrants  of  the  ^  plane  and 
the  Fourier  transforms  of  the  solutions  corresponding  to  the  complex  poles 
-  £   exist.   In  the  limit  as  the  loss  approaches  zero  (cr— v  0) ,  the  poles 
lie  on  the  real  axis.   The  surface  wave  solutions  still  satisfy  the 
homogeneous,  scalar  wave  equation  but  their  Fourier  transform  no  longer 
converges.   We  have  treated  the  limiting  case  where  evanishes  in  order  to 
be  able  to  solve  the  mode  equation.   This  approximation  is  justified  since, 
in  practice,  the  conductivity  o~  of  the  dielectric  is  extremely  small. 

We  know  from  the  solution  of  Maxwell  s  equations  for  the  source  free 
dielectric  rod  that  the  surface  wave  in  Region  I  is 

H^(p,z)  =  Ax  J^p)  e^°Z  =  Aj.  J]_(X1  g)  e^oz      (5. 26) 

Tangential  H  must  be  continuous  at  p  =  b  for  all  z,    so 


H^1  (b,z)  =  H^11  (b,z) 


Equating  (5.25)  and  (5.26)  at  p   =  b  yields 


1  ~  J1(x1)  An 


Therefore 


H  I/n   ^    a   V  }   T  (v     £>   i^oz  (5  =  27) 

H0  <p>z>  =  An  j^)  Ji(xi  V  e 
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We  obtain  EAp,z)    in  Regions  I  and  II  from  E  =  — - —  vy  H,  hence 
P  — iCoL 

?)H  I     £ 

V(p'z)  =  (I3l)  ^f  =  fe^0  V  (p'z)  (5-28) 

EpII(p,z)  =  (4^  V1  (p,z)      (529) 

The  field  components  (5.25),  (5,27),  and  (5.28),  and  (5.29)  will  be  used 
to  calculate  the  surface  wa.e  power  in  Chapter  6. 
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6.   THE  POWER  INTEGRALS  AND  EXCITATION  EFFICIENCY 

In  the  previous  section,  the  surface  wave  field  and  radiation 
field  generated  by  the  magnetic  current  source  were  determined.   It 
is  now  possible  to  calculate  the  corresponding  powers  and  obtain  the 
efficiency  of  the  source.   The  efficiency  with  which  the  source 
delivers  power  to  the  surface  wave  is  called  the  excitation  efficiency 
of  the  source.   Denoting  efficiency  by  the  symbol  T ,  it  is  defined  as 


S 

m  _  W_ 

;  -  T 
w 

s 

where  W  is  the  surface  wave  power 

T 
and   W  is  the  total  power  delivered  by  the  source „ 

12 
Goubau   has  proved  that,  for  lossless  surface  waveguides,  the 

radiation  and  surface  wave  fields  are  orthogonal  with  regard  to  power 

considerations.  In  this  case,  the  total  energy  delivered  by  the  source 

is  equal  to  the  sum  of  the  surface  wave  power  and  the  radiated  power. 

Since  we  are  considering  a  lossless  dielectric  rod,  the  orthogonality 

condition  holds  and  y  is  given  by 


ws  +  w* 

g 

where  W  is  the  surface  wave  power 

and   W  is  the  radiated  power. 

In  order  to  determine  T ,  we  proceed  with  the  calculation  of  W  and  W  . 


(6.1) 


(6.2) 
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6,1   Radiated  Power 

The  radiated  power  W  is  obtained  by  integrating  the  time  average 
Poynting  vector  over  the  surface  of  a  large  sphere  of  radius  r„ 

W  =J  jp"av  ■  da  =  j  |  Re  ¥  X  H*  •  da  (6.3) 


2  . 

where  da  =  r  cos  9  d  9  d<p. 

The  radiation  field,  which  has  components  E  and  H,,  is  independent  of 

0  and  is  symmetrical  about  the  plane  9=0,  so  (6,3)  reduces  to 

E 

2 
WR  =  2f      |      Re   E      H .*    (27Tr2   cos   9)    d9,  (6.4) 


Substitution  of    (5.19)    and    (5,20)    into    (6.4)   yields  the   integral    for 
the  radiated  power 


wR 


€        €      a   _ 

Id      {   b      } 

17 

-4-  r* 

(IT/2)    J 

o 

cos   9      I  F(9) I        d9„  (6,5) 


where 


J,  (k  aw) 

P(e)  = x    ° 


w  J   (k  bw)   H   (1)(k  b  cos  9)   -  £     cos   9   J   (k  bw)   H   (1)(k  b   cos  9) 
oolo  r  looo 


^j^T7 


2 

sin     9    „ 


A  discussion  of  the  integrand  of  (6.5)  is  worthwhile  at  this  point, 
We  note  that  the  magnitude  of  the  radiation  field  is  proportional  to 
|F(9)j,   As  9  approaches  UJ2,     |F(9) |  vanishes  so  the  radiation  pattern 


44 


has  a  null  along  the  axis  of  the  dielectric  rocL   The  influence  of  the 

source  dimension  (k  a)  on  the  radiation  field  and  W  is  contained  in 
o 

the  function  J,  (k  a  w)  which  is  the  numerator  of  F(0).   The  dielectric 
1   o 

rod  parameter  (k  b)  appears  only  in  the  denominator0   In  order  to  obtain 

R 

maximum  excitation  efficiency  for  any  k  b,  W  should  be  a  minimum,, 

W  will  be  a  minimum  if  J,  (k  a  w)  vanishes  or  is  quite  small  as  9 
1   o 

ranges  from  zero  to  7T/2;  this  occurs  when  J  (k  a  w)  passes  through  a 

I         2 
zero.   Consider  the  function  w  =\  €   -  sin  0  with  €   =  2t56.   As  9 

N   r  r 

takes  on  values  0  ^  Q  ^  ^12,    w  has  the  range  1„25  ^  w  £   1.6.   Assume 

an  average  value  for  w,  w   =  1.42.   We  should  expect  the  most  efficient 

source  to  have  the  dimension  k  a«^  3.83/w   =  2.7,  where  3.83  is  the 

o        '  av 

first  zero  of  the  J  function.   This  analysis  is  verified  in  the 
efficiency  data  presented  in  the  latter  part  of  this  section. 

The  integral  portion  of  (6.5)  was  evaluated  numerically  on  the 
University  of  Illinois  digital  computer.   Since  the  denominator  of 
F(9)  is  unbounded  as  9  approaches  1FJ2,    the  upper  limit  of  integration 
was  selected  as  88  and  the  range  of  integration  was  subdivided  into 
44  equal  intervals  of  2  each.   The  computer  integration  routine 
evaluated  the  integral  by  Simpson's  rule.   The  computer  program  is 
discussed  more  fully  in  Appendix  Co   The  constant 


€    £  a  2 

\Jl  K    b  ' 


g 
which  multiplies  the  integral  in  (6.5)  also  occurs  in  W  and  may  be 

cancelled  in  the  evaluation  of  T .   The  quantity  that  was  determined 

on  the  computer  and  which  was  used  to  calculate  T  is  the  normalized 

integral 

88  £ 
90  2 

2 


nR  =  571   P   cos  e  I  F(e)  i  de  (66) 
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where 


wR 


€   €  a  _ 

a  v  b  ' 

o 


p 
N  was  computed  as  k  a  varied  in  discrete  steps  of  0.2  over  the  range 

0  C   k  a  ^  k  b  for  the  selected  value  of  k  b  under  consideration.   The 
o  -  o  o 

S 
results  are  given  in  Table  II,  along  with  N  and  the  excitation 

efficiency  **p  . 

6.2  Surface  Wave  Power 

The  total  power  in  the  surface  wave  launched  by  the  magnetic  current 

source  is  determined  in  the  following  manner.   Consider  an  infinite 

transverse  plane  normal  to  the  z  axis  located  at  some  position  z 

positive.   The  power  propagating  in  the  positive  z  direction  is  obtained 

by  integrating  the  time  average  Poynting  vector  over  the  surface  of  the 

infinite  plane.   Because  of  symmetry,  the  power  propagating  in  the 

negative  z  direction  is  equal  to  the  power  propagating  in  the  positive 

z  direction.   Therefore,  the  total  surface  wave  power  is  equal  to  twice 

the  power  in  the  positive  z  direction.   Since  the  field  components 

E  ,  H.  are  different  for  Regions  I  and  II,  separate  integrations  are 

necessary  for  the  regions  0  C  p  c   b  and  p  y   b.   Consider  the  positive 

z  direction.   We  denote  the  power  over  the  cross  sectional  area  of  the 

rod  as  W  and  W   applies  to  region  II,  which  extends  to  infinity.   The 

a 

total  surface  wave  power  W  is  given  by 


WS  i=  2  (W1  +  W   ).  (6.7) 


Consider  region  I  where  0  ^  p  £_  b,    the   integral   for  W     is 


27T     p=b 

rs 

o     p=o 


W1   =       j       J  1/2  Re  E    *   H,1*  da  (6.8) 
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where  da  =  pdpdcf). 

Substituting  (5.38)  and  (5.27)  into  (6,8),  we  obtain 


w  ■  ^    tAn  j^y]      J    p  h2  (\  f)  <p.  (6.9) 

p=o 

The   value  of  the  definite  integral  in  (6.9)  is  given  by  McLachlan,13 
so  (6.9)  reduces  to 

wI  =  r  fer>  [An  j^yll    V<  Jo2(xi)  +  Ji2  (xi>  -  ^  VV  VxiY   («.io) 

tor  region  II  where  p   >  b,  W11  is  given  by 


27T       p=o> 

w"=/  /  i'2  ReEPn  VI%ai  (6-u> 

o       p=b 
Substituting    (5.29)    and    (5.25)    into    (6,11)   yields 

p=oo 


fi--&*ii  r  p»ka'*&* 


p=b 


(6.12) 


he  definite  integral   in   (6.12)    is  also  given  by  McLachlan,14 


btain  W11  as 


w"  ■  <£f>  A„2  4U02  <6>  -  kx2  (I)  +  I  ko  (?)  Kl  <e>j 

le  constant  A^  is  defined  in  Eq.  (5.25),   We  may  now  substitute 
5»10),  (6.13),  and  the  value  of  A   into  (6.7)  to  obtain  WS.   After 


(6.13) 


considerable  manipulation  of  constants  one  obtains  the  expression 
for  the  total  surface  wave  power  which  is 


47 


where 


n  x  \ 

SI    o 


o 


K,  (I) 


Nr<ry]    &«  (V+Ji  <V  -  r  VV  VV> 


r   1 x  1 


(6.14) 


"(4)2+1l  r     a    X!  n 

3— ■jJ0(X1)K1(6)  +  [€r(^)  +  (^Jo(X1)Ko(e)  + 


[Ko2(e)-Kl2(e)  +  f  ko<£)  ^«a/ 


(€r  -  1)  ji   (x1)  1^  (I) 


2   V 

Ji  <Fb  V 

o 


Inspection  of  the  expression  for  N  shows  that  the  source  parameter 

k 


(k  a)  appears  only  in  the  argument  of  J. 


i2(k\X,) 


The  remaining 


o 

portion,  which  is  rather  formidable,  is  a  function  6  and  the  mode 

r 

solution  (£,X  ) „   It  has  a  constant  value  for  any  €   and  k  b  which, 

t  S 

of  course,  yields  the  eigenvalues  5  and  X  .   N  is  a  function  of  the 


source  dimension  only  by  the  term  J,   1  - — -  X, 
J  y  1   Vk  b  V 

o 


It  should  be  noted, 


however,  that  a/ b  appears  in  the  constant  which  multiplies  N  to 

S  S 

give  W  .   Selecting  k  b  with  €   =  2.56,  N  was  calculated  as  a 
o        r 

function  of  k  a  using  values  of  X  !\  „  £,  and  Xn  from  Table  I 
o  gl  o'         1 

computed  values  of  N  are  given  in  Table  II . 


The 
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6.3  Excitation  Efficiency 

The  excitation  efficiency  rp  for  the  magnetic  current  ring  on  a 
lossless  dielectric  rod  is  given  by 


w  +  w 


£    £ra  2 


When  (6.6)  and  (6.14)  are  substituted  into  (6.2)  the  constant 

\ 

may  be  cancelled,  so  that  (6.2)  reduces  to 

NS 
NS  +  N* 

Excitation  efficiency  was  calculated  according  to  (6.15)  using  values 

S       R 
of  N  and  N  obtained  for  the  particular  source  dimension  k  a  and  rod 

o 

S   R 
parameter  k  b.   Table  II  presents  N  ,  N  ,  and  ^  for  each  k  a  and  k  b„ 

Figures  8  and  9  show  curves  of  efficiency  as  a  function  of  k  a  for 

o 

the  selected  values  of  k  b.   The  curves  are  graphs  of  the  data  presented 

in  Table  II.   The  analysis  indicates  that  a  magentic  current  ring  of 

the  proper  size  will  provide  very  efficient  excitation  of  the  E   mode. 

Efficiencies  greater  than  90%  are  predicted  for  k  a  approximately  equal 

o 

to  2.6,  which  agrees  very  well  with  the  rough  estimate  of  optimum 
dimensions  that  was  made  in  Section  6.1.   Experimental  measurements  of 
efficiency  were  performed  in  the  laboratory  as  a  check  on  the  foregoing 
analysis.   The  experimental  work  is  described  in  Chapter  7. 


TABLE    II 


49 


k   b 
o 


k   a 


NR 


2.2 


2.6 


3.0 


.2 

,0179 

.0673 

.2101 

.4 

.0684 

.2,507 

.2144 

.6 

.1427 

.5010 

.2217 

.8 

.2278 

.7517 

.2326 

1.0 

.3087 

.9379 

.2476 

1.2 

.3719 

1.0130 

.2685 

1.4 

.4066 

.9621 

.2971 

1.6 

.4077 

.8026 

.3369 

1.8 

.3755 

.5784 

.3936 

2.0 

.3163 

.3469 

.4769 

2.2 

.2403 

.1604 

.5997 

.2 

.0320 

.1162 

.2159 

.4 

.1228 

.4337 

.2207 

.6 

.2582 

.8659 

.2297 

.8 

.4171 

1 . 2969 

.2433 

1.0 

.5744 

1.6129 

.2626 

1.2 

.7064 

1.7357 

.2892 

1.4 

.7928 

1.6380 

.3261 

1.6 

.8214 

1.3540 

.3776 

1.8 

.7891 

.9625 

.4505 

2.0 

.7020 

.5620 

.5554 

2.2 

.5741 

.2448 

.7011 

2.4 

.4249 

.0692 

.8599 

2.6 

.2764 

.0481 

.8518 

.2 

.0400 

.1169 

.2549 

.4 

.1544 

.4369 

.2611 

.6 

.3278 

.8753 

.2725 

.8 

.5370 

1.3165 

.2879 

1.0 

.7536 

1.6489 

.3137 

1.2 

.9493 

1.7897 

.3466 

1.4 

1.0984 

1.7099 

.3911 

1.6 

1.1817 

1.4377 

.4511 

1.8 

1.1902 

1 . 0474 

.5319 

2.0 

1,1234 

.6371 

.6381 

2.2 

.9920 

.2988 

.7685 

2.4 

.8136 

.0962 

.8943 

2.6 

.6115 

.0481 

.9271 

2.8 

.4109 

.1302 

.7594 

3.0 

.2349 

.2859 

.4510 

TABLE    II    (CONTINUED) 
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k  b 
o 


k  a 


NR 


3.4 


3.8 


.2 

.0371 

.1189 

.2378 

.4 

.1442 

.4443 

.2450 

.6 

.3090 

.8917 

.2573 

.8 

.5130 

1 . 3470 

.2758 

1.0 

.7329 

1.6951 

.3019 

1,2 

.9438 

1.8523 

.3375 

1.4 

1.1228 

1.7866 

.3859 

1.6 

1 . 2501 

1.5238 

.4507 

1.8 

1.3124 

1.1342 

.5364 

2.0 

1.3038 

.7142 

.6461 

2.2 

1.2271 

.3575 

.7743 

2.4 

1.0913 

.1306 

.8931 

2.6 

.9076 

.0583 

.9396 

2.8 

.7087 

.1212 

.8540 

3.0 

.5032 

.2667 

.6536 

3.2 

.3164 

.4290 

.4244 

3.4 

.1646 

.5487 

.2308 

.2 

.0334 

.1561 

.1763 

.4 

.1304 

.5843 

.1825 

.6 

.2816 

1.1733 

.1936 

.8 

.4723 

1.7725 

.2104 

1.0 

.6845 

2.2325 

.2347 

1,2 

.8984 

2.4437 

.2688 

1.4 

1.0933 

2.3639 

.3162 

1.6 

1.2511 

2.0228 

.3821 

1.8 

1.3583 

1.5159 

.4726 

2.0 

1 . 4048 

.9668 

.5923 

2.2 

1.3873 

.4983 

.7357 

2.4 

1.3088 

.1967 

.8693 

2.6 

1.1772 

.0946 

.9255 

2.8 

1.0056 

.1674 

.8573 

3.0 

.8096 

.3485 

.6991 

3.2 

.6073 

.5530 

.5234 

3.4 

.4154 

.7052 

.3767 

3.6 

.2492 

.7584 

.2473 

3.8 

.1209 

.7048 

.1464 
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TABLE  II  (CONTINUED) 


k  b 
o 

4  .  2  .  0967 

1010 
1085 
1202 
1378 
1632 
2009 
2569 
3414 
4706 
6389 
8079 
8763 
7861 
6232 
4748 
3646 
,2859 
,2273 
,1771 
,1237 


k  a 
o 

NS 

NR 

.2 

.0290 

.2710 

.4 

.1138 

1.0130 

.6 

.2471 

2.0307 

,8 

.4184 

3.0632 

1.0 

.6140 

3.8423 

1.2 

.8177 

4.1927 

1.4 

1.0133 

4.0300 

1.6 

1.1854 

3.4293 

1.8 

1.3205 

2.5469 

2.0 

1.4326 

1.6114 

2.2 

1.4671 

.8292 

2.4 

1.4211 

.3379 

2.6 

1.3465 

.1901 

2.8 

1.2259 

.3336 

3.0 

1.0687 

.6461 

3.2 

.8876 

.  9817 

3.4 

.6967 

1 . 2140 

3.6 

.5096 

1.2727 

3.8 

.3396 

1.1546 

4.0 

.1973 

.9168 

4.2 

.0912 

.6461 
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7.   EXPERIMENTAL  INVESTIGATION 

7.1   The  Source  Representation  and  Measurement  Method 

The  magnetic  current  ring  was  purposely  chosen  for  the  theoretical 
problem  because  it  is  a  source  which  has  a  simple  physical  realization. 
Consider  Fig.  1.   By  symmetry  the  magnetic  current  filament  generates  a 
field  which  has  no  electric  component  in  the  plane  z  =  0.   In  the 
plane  of  the  source  the  only  non-zero  component  of  E  is  E  ,  which  is 
normal  to  the  plane.   Since  tangential  E  vanishes  over  the  plane,  an 
infinitely  large  metal  sheet  may  be  placed  at  the  source  position 
without  disturbing  the  field.   One  may  simulate  the  source  by  a  very 
narrow  annular  slot  in  a  large  conducting  sheet  or  ground  plane.   Then, 
of  course,  the  dielectric  rod  extends  from  the  ground  plane  in  only  one 
direction.   The  annular  slot  may  be  properly  illuminated  by  a  circular 
or  coaxial  waveguide.   If  the  rod  is  terminated  in  a  resistance  card 
load  which  results  in  negligible  reflection  of  the  surface  wave,  then 
the  finite  rod  excited  by  an  annular  slot  in  a  ground  plane  simulates  in 
a  half  space  the  theoretical  problem  of  Fig.  1.   A  photograph  of  the 
experimental  structure  assembled  in  the  laboratory  is  given  in  Fig.  10. 

In  order  for  the  slot  to  be  a  good  approximation  of  the  filament, 
the  slot  width  must  be  small  compared  to  the  wavelength  and  the 
electric  field  in  the  slot  must  be  radial  and  uniform  about  the 
circumference.   It  is  assumed  that  the  mean  radius  of  the  slot  corresponds 
to  the  radius  of  the  infinitesimal  current  filament.   The  slot  may  be 
illuminated  from  a  circular  waveguide  propagating  the  TM   mode,  or 
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FIGURE  10.   VIEW  OF  THE  DIELECTRIC  ROD  AND  GROUND  PLANE 
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from  a  coaxial  line  excited  in  the  usual  TEM  mode.  The  coaxial  line 
is  preferable  since  the  center  conductor  provides  a  convenient  means 
of  supporting  the  circular  disc  which  is  the  ground  plane  within  the 
annular  slot. 

The  efficiency  with  which  the  slot  excites  the  surface  wave  may  be 

5,6 
measured  using  Deschampc  ,.iethod    for  determining  the  scattering 

matrix  coefficients  of  a  waveguide  junction.   The  transition  from  the 

coaxial  line  to  the  dielectric  rod  waveguide  is  considered  as  a  two 

port  waveguide  junction.   Assuming  that  the  dielectric  rod  and  coaxial 

line  are  lossless,  then  the  dissipative  attenuation  of  the  junction 

results  entirely  from  the  power  radiated  into  space  by  the  slot.   If 

all  of  the  power  delivered  to  the  slot  is  converted  to  the  surface 

wave  mode,  then  there  is  no  power  lost  as  stray  radiation  and  the 

junction  is  100%  efficient.   This  method  has  been  used  previously  to 

measure  the  launching  efficiency  of  sources  placed  on  a  dielectric 

15 
image  line,, 

The  experimental  procedure  for  measuring  the  efficiency  of  the 

junction  is  illustrated  in  Fig.  11. 


(Dete^Vor  J 
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Coaxial    <5Jo4+ed 


Junction 
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FIGURE   11.       REPRESENTATION  OF  A  TWO   PORT   JUNCTION 
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The  short  circuit  plate  which  terminates  the  dielectric  rod  is 
placed  at  four  arbitrary  but  equally  spaced  positions.   The  reflection 
coefficient  corresponding  to  each  position  of  the  short  is  measured 
in  the  coaxial  line  and  plotted  on  the  reflection  coefficient  plane „ 
The  measured  points  define  a  circle  which  is  the  image  of  the  unit 
circle  in  the  output  guide .   Using  the  graphical  constructions  derived 
by  Qeschamps,  one  may  determine  the  crossover  point,  inconocenter ,  and, 
ultimately,  the  scattering  matrix  coefficients  which  characterize  the 
junction. 

The  graphical  construction  used  to  determine  the  iconocenter  may 
be  avoided  if  a  matched  load  is  available  for  the  output  guide,  since 
the  coefficient  S   and  the  iconocenter  are  the  same  point  on  the 
reflection  coefficient  plane.   If  one  terminates  the  output  guide  in 
a  matched  load,  the  reflection  coefficient  measured  in  the  input  guide 
is  S   and,  therefore,  the  iconocenter. 

The  so-called  vVmatched"  efficiency,  no  ,  is  defined  as 

Pm  !S12|2 

P      ~  I  1 2  * 

1   l~ls»l 

where  P  is  the  power  delivered  to  a  matched  load  (infinitely  long  guide) 
m 

and  P.  is  the  power  delivered  to  the  junction  when  the  output  guide 
is  terminated  by  a  matched  load.   We  see  that  T   is  identical  to  the 
excitation  efficiency  defined  by  Eq.  (6.1).   The  relation  between  f 
and  the  scattering  coefficients  is  given  in  Eq.  (7,1).   In  order  to 
calculate  T  one  needs  only  the  iconocenter  and  the  locus  of  the  image 
circle.   A  single  measurement  for  the  matched  load  termination  determines 
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the  iconocenter,  and  three  measurements  for  arbitrary  positions  of 
the  short  circuit  plate  define  the  image  circle,   Mathis   gives  a 
convenient  formula  f or  HP  which  is  equivalent  to  the  relation  (7.1). 
His  formula  is 

2D1D2 
T-  5-  (7.2) 

<VV  (1  "  D3  > 

where  D  ,  D  ,  and  D  are  the  lengths  defined  in  Fig.  12.   The  origin 

of  the  reflection  coefficient  plane  is  the  point  0,  the  iconocenter 

is  0   ,  and  the  center  of  the  image  circle  is  C.   Formula  (7.2)  was 

used  to  calculate  the  efficiency  for  all  of  the  experimental  measurements. 

Before  the  results  are  presented,  a  brief  description  of  the  experimental 

apparatus  and  coaxial  exciter  will  be  given. 

7.2  Experimental  Equipment 

A  two  inch  diameter  polystyrene  rod  was  selected  for  the  experimental 
work.   This  choice  followed  from  two  main  considerations.   First,  a  two 
inch  rod  can  stand  on  end  without  additional  support.   This  method  of 
vertical  mounting  avoids  the  usual  problem  of  how  to  support  a  surface 
wave  structure  without  disturbing  the  field.   Second ,  it  is  easier 
to  fabricate  annular  slots  for  a  large  rod  than  for  a  small  rod. 
Table  III  lists  the  wavelength  X  ,  frequency  f ,  the  ratio  X  (X  ,  and 
X  for  a  two  inch  polystyrene  rod,  assuming  €   =2.56. 

Using  Deschamps7  method,  one  obtains  the  efficiency  of  the  entire 
transition  between  the  input  reference  plane  and  the  output  reference 
plane,  that  is,  between  the  measuring  probe  and  the  short  circuit 
termination.   Since  we  are  interested  in  measuring  the  dissipative 
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Unit  Circle 


FIGURE  12 .   PLOT  OF  THE  IMAGE  CIRCLE  ON  THE  REFLECTION 
COEFFICIENT  PLANE 


TABLE  III 

k  b 
o 

X 

o 

f 

g/  o 

X 

g 

centimeters 

megacycles 
sec 

centimeters 

2  2 

7  25 

4132 

0.969 

7.03 

2,6 

6.14 

4884 

0.904 

5.55 

3.0 

5.32 

5635 

0.840 

4.47 

3.4 

4.69 

6387 

0.791 

3.71 

3.8 

4.20 

7138 

0.756 

3.17 

4.2 

3.80 

7889 

0.730 

2.78 

attenuation  of  just  the  annular  slot,  it  is  essential  that  the  rest  of 
the  transition  shall  introduce  only  negligible  attenuation.   For  the 
present  we  shall  ignore  the  dielectric  rod  loss  since  it  is  very  small. 
It  follows,  then,  that  a  low  loss  exciter  must  be  constructed  to 
illuminate  the  annular  slot.   The  annular  slot  presents  a  very  low 
conductance  and  capacitive  susceptance  to  the  exciting  waveguide. 
Therefore,  the  feed  waveguide  must  be  a  low  impedance  line  with  some 
means  of  tuning  out  the  capacitive  susceptance  of  the  slot  at  the 
ground  plane  position.   Figure  13  shows  a  cross  sectional  view  of  the 
low  impedance  coaxial  line  which  was  constructed  for  this  purpose. 
The  inner  diameter  of  the  outer  wall  of  the  coax  was  17/8  inches 
and  the  inner  conductor  diameter  was  1  l/4  inches,  which  yields  a 
24  ohm  line.   A  two  wavelength  tapered  section  transformed  the  24  ohm 
line  to  a  standard  50  ohm,  type  N  connector.   Two  polystyrene  rings 
centered  the  inner  conductor  within  the  cylinder.   The  ring  located 
at  the  tapered  section  was  one-half  wavelength  long,  so  that  no  impedance 
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discontinuity  would  result.   The  ring  at  the  ground  plane  position 

was  a  quarter-wavelength  transformer  which  would  match  the  24  ohm  line 

to  a  9.5  ohm  resistance  load.   A  circular  tuning  disc  or  washer  was 

placed  on  the  end  of  the  center  conductor  to  provide  a  series  inductive 

reactance  to  cancel  the  capacitive  susceptance  of  the  slot.   The  circular 

disc  which  formed  the  inner  boundary  of  the  annular  slot  was  fixed  to 

the  tuning  washer  and  the  coax  center  conductor  by  a  special  mounting 

screw.,   Since  the  neighborhood  of  the  annular  slot  and  tuning  washer  is 

a  resonant  cavity,  all  parts  were  silver-plated  to  minimize  losses,. 

A  family  of  discs  and  rings  were  fabricated  to  permit  varying  the 

annular  slot  radius  from  1/2  inch  to  7/8  inch  while  maintaining  the 

slot  width  constant  at  1/8  inch.   The  1/8  inch  slot  width  corresponds 

to  .067  X  and  .076  X  at  k  b  equal  to  3.4  and  3.8,  respectively, 
o  o     o 

Recalling  that  the  radius  of  the  dielectric  rod  is  one  inch,  it  is 

convenient  to  express  the  slot  radius  in  the  normalized  from  a/b; 

then  the  source  dimension  k  a  for  any  k  b  is  given  by  k  b  (a/b)„   Six 

o  o  o     ' 

slots  were  constructed  for  the  measurements.   Table  IV  gives  the 

normalized  slot  radius  a/b,  and  the  corresponding  dimension  k  a  at 

the  two  values  of  k  b  which  were  used  for  the  measurements, 
o 

A  slot  1/8  inch  wide  and  two  inches  long  was  milled  in  the  side 
wall  of  the  exciter  so  that  the  standing  wave  ratio  could  be  measured. 
After  a  particular  slot  was  mounted  on  the  exciter,  various  tuning 
washers  were  tested  until  one  was  found  which  reduced  the  VSWR  in  the 
exciter  to  less  than  two  when  the  dielectric  rod  was  terminated  with 
a  matched  load.   The  image  circle  was  determined  from  measurements 
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TABLE  IV 


Normalized  Slot  o  o 

Radius  a/b  when  k  b=3.4  when  k  b^3.8 

o  o 


0.50  1.70  1.90 

0,625  2.12  2.38 

0.687  2.34  2.61 

0.75  2.55  2.85 

0.812  2.76  3.08 

0.875  2.98  3.32 


made  in  a  coaxial  slotted  line  which  was  connected  directly  to  the 

exciter  through  a  type  N  elbow  connector.   A  view  of  the  slotted  line, 

coaxial  exciter ,  and  the  underside  of  the  ground  plane  is  given  in 

Fig.  14.   A  block  diagram  of  the  waveguide  apparatus  used  for  the 

measurements  is  given  in  Fig.  15.   When  Desehamps'  procedure  was 

carried  out,  the  dielectric  rod  was  terminated  with  a  six  inch 

diameter  plate,  machined  from  1/16  inch  brass  and  silver  plated.   The  ground 

plane,  which  may  be  seen  in  Fig.  10,  was  60  inches  square.   For  most 

of  the  measurements  the  dielectric  rod  was  40  centimeters  long, 

although  in  some  instances,  the  length  was  increased  to  134  centimeters. 

Finally,  it  should  be  mentioned  that  an  external  probe  was  constructed 

to  measure  the  standing  wave  ratio  of  the  field  along  the  dielectric 

rod.   It  is  shown  schematically  in  Fig.  15.   It  was  used  to  measure 

the  guide  wavelength,  X  „  when  the  rod  was  operated  at  k  b  equal 

to  3.4  and  3.8.   In  each  case,  the  measured  guide  wavelength  corresponded 

exactly  to  the  value  predicted  in  Table  III.   A  matched  load  was  made 

from  four  tapered  lengths  of  300  ohm  resistance  cord  which  were 
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mounted  symmetrically  on  the  rod.   The  matched  load  is  visible  in 
Fig  10   The  VSWR  on  the  dielectric  rod  was  less  than  1.1  when  the 
rod  was  terminated  with  this  load, 
7  3   Experimental  Results 

When  efficiency  was  measured  as  shown  in  Fig.  11,  the  result  obtained 
was  not  precisely  the  excitation  efficiency  of  the  annular  slot;  instead 
it  was  the  efficiency  of  the  entire  transition  between  the  measuring 
probe  and  the  short  circuit  plate.   Although  the  losses  in  the  system 
were  small,  the  measured  efficiency  was  reduced  somewhat  by  dielectric 
rod  attenuation  and  by  the  loss  between  the  measuring  probe  and  the 
ground  plane.   We  shall  discuss  these  two  effects  separately. 

First,  consider  the  dielectric  rod  attenuation  which  caused  the 
two  port  network  to  be  terminated  in  a  lossy  short  circuit.   The  lossy 
rod  acts  as  a  reflection  coefficient  transformer,  so  that  the  magnitude 
of  the  reflection  coefficient  at  the  ground  plane,  IfL  is  given  by 


h 


?   "''"  (7.3) 


where  a  is  the  attenuation  constant  of  the  wave  on  the  dielectric 
rod  and  z  is  the  length  of  rod  between  the  ground  plane  and  the  short 

ircuit  plate.   It  is  assumed,  of  course,  that  a  perfect  short  terminates 
the  rod  yielding  a  reflection  coefficient  of  unity  at  that  position.   We 
wish  to  determine  the  efficiency  of  the  slot  while  excluding  the  effect 
of  the  lossy  rod.   Consider  Fig.  16(a),  which  shows  [s']  ,  the  four 
Pole  to  be  measured;  it  is  located  between  the  slotted  line  and  the 
lossy  short  circuit.   Figure  16(b)   is  an  alternate  representation  of 
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Lossy      Short 


51.      %* 


I- 


(a) 


Pev-fec.+  Shovl- 


(b) 


FIGURE  16.   MEASUREMENT  OF  A  JUNCTION  WITH  A  LOSSY  SHORT  CIRCUIT 


Fig.  16(a)  in  which  the  reflection  coefficient  transformer  I  P|has 
been  associated  with  [s» ]  to  form  the  fictitious  four  pole  [s]  which 
is  terminated  with  a  perfect  short,   when  one  measures  the  four  pole, 
assuming  a  perfect  short,  the  composite  four  pole  [s]  is  obtained. 
Assuming  that  |  P  |is  known,  the  parameters  of  [S»  ]  are  given  by] 


17 


S   » 
11 


11 


S   ' 
22 


_22 


12 


"12 


Referring  to  Eq.  (7.2),  it  is  evident  that  the  efficienty,^ f  ,  of 
[S>  ]  is  related  to  T  of  [S]  by 


T«  = 


T 


(7.4) 


(7.5) 


is  the  slot  efficiency  which  would  be  measured  if  the  dielectric 
rod  were  lossless,  and  Eq.  (7.5)  gives  the  relation  between  the  measured 


6!' 


p 


efficiency  T»  and  T  '  . 

In  order  to  apply  (7.5)  to  the  experimental  measurements,!1   I  was 
obtained  by  measuring  the  efficiency  of  the  same  source  with  two 
different  lengths  of  dielectric  rod.   The  distance  from  the  ground 
plane  to  the  short  circuit  plate  was  40  centimeters  and  134  centimeters 
for  the  two  measurements.   Referring  to  Fig.  17,  let T  and  f   represent 
this  efficiencies  which  would  be  measured  when  the  short  was  located 
at  the  specified  positions. 


GjirOund  Plane. 


Die-le+vic  Rod  — 7 
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-.94n 
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FIGURE  17,   ARRANGEMENT  FOR  MEASURING  ROD  ATTENUATION 


From  (7,3)  and  (7.5)  it  follows  that 


£-•*»  (7.6) 

(o 


where  z  is  equal  to  0,94  meters,, 

Substituting  the  measured  values  of  T  and  HP  into  (7.6),  one 
obtains  the  value  of  the  attenuation  constant  a.   Having  determined  a, 
it  follows  from  (7,5)  that 

T»  =  Te2aZ  (7.7) 

where  z  is  now  0,40  meters,  '"P '  may  be  thought  of  as  the  measured 

efficiency  referred   to  the  ground  plane  position.   Equation  (7,7) 

defines  the  adjustment  which  is  applied  to  the  measured  efficiency,  T, 

to  account  for  the  dielectric  rod  attenuation.   The  value  of  a  was 

measured  when  k  b  was  equal  to  3.4  and  3.8,   The  results  are 
o 


-    ,0692, 


k  b  =  3o4 


,0716 


k  b  =  3.8 


(7.8) 


o  o 

Before  (7.7)  is  applied  to  the  experimental  measurements,  the  attenuation 

between  the  measuring  probe  and  the  ground  plane  will  be  discussed. 

The  attenuation  of  the  transition  between  the  measuring  probe 

and  the  ground  plane  results  from  losses  in  the  type  N  connector 

and  the  coaxial  exciter.   The  loss  is,  of  course,  a  function  of  the 

standing  wave  ratio  in  the  transition.   When  a  particular  slot  is  being 

measured ,  the  VSWR  changes  with  the  short  circuit  position  because  the 

junction  is  lossy.   Therefore,  the  transition  loss  is  different  for 
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each  position  of  the  short,  and  it  is  impracticable  to  correct  every 

standing  wave  measurement  performed.   One  may,  however,  establish  the 

maximum  loss,  or,  correspondingly,,  the  minimum  efficiency  of  the 

transition.   This  was  accomplished  by  mounting  a  two  inch  diameter, 

circular  waveguide  perpendicular  to  the  ground  plane  and  concentric  to 

the  annular  slot,  and  measuring  the  efficiency  by  Deschamps'  method. 

The  circular  guide  was  fitted  with  a  movable  shorting  plunger.   The 

annular  slot  excites  the  TM __  mode  in  the  guide.   The  schematic 

representation  of  the  two  port  junction  is  shown  in  Fig.  18. 


I  Detector  j 


Qround    Plane 


Circular 
Waveguide 


FIGURE  18.      ARRANGEMENT   FOR  MEASURING  THE  TRANS I ST ION 
EFFICIENCY 
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Assuming  that  the  circular  guide  was  lossless,  the  measured 
efficiency,  T,  is  the  efficiency  of  the  transition  from  the  slotted 
line  to  the  annular  slot.   It  is  now  possible  to  adjust  the  measured 
slot  efficiency  to  account  for  both  the  transition  loss  and  the 
dielectric  rod  attenuation.   Let  Tq   denote  the  corrected  slot  efficiency; 
it  is  given  by 

2az 


T 


TL_  Tel 

TT~   TT 


(7.9) 


where  Tis  the  efficiency  measured  with  the  lossy  rod  present,  and 
TT  is  the  transition  efficiency.   It  should  be  noted  that  T  provides 
a  maximum  correction  to  T  since,  in  most  cases,  the  transition  efficiency 
will  be  greater  than  the  value  Xp-   This  results  because  standing  wave 
ratios  are  higher  when  "^   is  measured  than  when  T  is  determined.  T 
was  measured  when  kQb  was  equal  to  3.4  and  3.8.   The  results  are 


T 


0.939, 


T 


kb=3,4 
o 


=  0.934 


(7.10) 


b  =  3.8 


Substituting   (7.8)   and    (7.10)    into    (7.9)   yields 


*   1.125  T,      T 


k  b  =  3.4 
o 


=   1.132T 


k  b  =  3.8 
o 


(7.11) 


The  measured  efficiencies  were  corrected  according  to  Eqs.  (7.11). 
tote  that  the  correction  factor  is  approximate.   It  was  obtained 
Ln  order  to  fix  the  order  of  magnitude  of  the  system  losses.   One 
lay  see,  however,  that  the  excitation  efficiency  of  a  slot  is  somewhat 
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greater  than  the  measured  value. 

Six  slots  were  measured  in  the  laboratory;  the  results  are  presented 
in  Table  V,  which  includes  the  measured  efficiency  °f  and  the  corrected 
efficiency  f   for  each  source  k  a„   Measurements  were  performed  at 
frequencies  of  6387  megacycles/sec  and  7138  megacycles/sec  which 
result  in  k  b  equal  to  3.4  and  3.8,  respectively.   The  data  of  Table  V 
have  been  plotted  in  Figs.  19  and  20  for  comparison  with  the  theoretical 
curves  of  efficiency.   The  very  close  agreement  between  the  experimental 
points  and  the  theoretical  curve  is  evident.   In  all  cases  the  measured 

TABLE  V 

k  b  k  a  Measured  Corrected 

Efficiency  Efficiency 

T  Tc 

3  4  1.70  0.44  0.495 

2.12  0.65  0.73 

2.34  0.76  0.855 

2.55  0.835  0.94 

2.76  0.77  0.867 

2.98  0.65  0.732 

3  8  1.90  0.48  0.544 

2.38  0.775  0.88 

2.61  0.85  0.96 

2.85  0.77  00872 

3.08  0,63  0.714 

f  3.32  0.44  0.498 

efficiency,  which  included  the  system  losses,  was  within  10  percent  of 

the  efficiency  predicted  by  theory.   The  experimental  measurements 

verify  that  an  excitation  efficiency  of  approximately  95  percent  may 

be  obtained  from  an  annular  slot  of  dimension  k  a  =  2.6. 
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7.4  Elimination  of  the  Ground  Plane 

The  large  ground  plane  was  necessary  in  order  for  the  annular  slot 

to  properly  represent  the  theoretical  source.   A  point  of  practical 

interest,  however,  is  the  efficiency  of  the  annular  slot  in  the  presence 

of  a  small  ground  plane.   This  effect  was  investigated  by  measuring 

slot  efficiencies  as  the  size  of  the  ground  plane  was  reduced.   Surprisingly, 

the  ground  plane  had  little  effect  on  the  efficiency.   The  ground  plane 

was  reduced  to  a  10  inch  diameter  flange  and  no  appreciable  change  in 

efficiency  was  noted.   Consequently,  the  ground  plane  was  eliminated 

entirely  so  that  the  experimental  structure  consisted  of  the  dielectric 

rod  mounted  on  the  coaxial  exciter.   The  efficiency  of  each  slot  was 

measured  at  k  b  equal  to  3.4  and  3.8.   The  results  are  presented  in 
o 

Table  VI  which  indues,  for  comparison,  the  efficiency  that  was  measured 
with  the  large  ground  plane.   It  may  be  seen  that  elimination  of  the 
ground  plane  had  little  effect  on  the  launching  efficiency  of  this 
source.   The  annular  slot  in  the  end  of  the  coaxial  section  is  a  very 
efficient  exciter  of  the  E   mode. 


k  b 
o 


3,4 


3.8 


TABLE  IV 

k  b 

Efficiency  T 

Efficiency  ^ 

o 

with  Ground 

Plane 

No  Ground  Plane 

1,70 

0.44 

0.428 

2.12 

0.65 

0.647 

2.34 

0.76 

0.779 

2.55 

0.835 

0.847 

2.76 

0.77 

0.829 

2.98 

0.65 

0.696 

1.90 

0.48 

0.545 

2.38 

0.775 

0.776 

2.61 

0.85 

0.835 

2.85 

0.77 

0.776 

3.08 

0.63 

0.63 

3.32 

0.44 

0.413 

76 
CONCLUSIONS 

The  theoretical  analysis  predicted  that  a  circular  filament  0.827 
wavelengths  in  diameter  would  launch  the  E   mode  with  an  efficiency 
of  approximately  95  percent.   The  analysis  also  showed  that  excitation 
efficiency  is  a  critical  function  of  the  source  dimension,  and  high 
efficiency  is  obtained  only  for  a  limited  range  of  source  diameters. 
It  is  evident  that  a  filament  which  is  greater  in  diameter  than  the 
dielectric  rod  is  a  poor  exciter  of  the  surface  wave. 

The  uniformly  illuminated  annular  slot  was  an  excellent  physical 
representation  of  the  theoretical  source.   The  measured  slot  launching 
efficiencies  were  in  close  agreement  with  the  results  predicted  by 
the  theory.   It  was  determined  that  the  annular  slot  efficiency  is 
essentially  independent  of  the  size  of  the  ground  plane.   Very  good 
efficiency  may  be  obtained  using  a  simple  launcher  such  as  the  coaxial 
exciter. 

In  order  to  realize  the  high  efficiency  of  the  annular  slot,  it 
is  necessary  to  illuminate  the  slot  from  a  good  ^nwtched"  transition. 
In  the  laboratory,,  the  predicted  source  efficiencies  were  not  measured 
until  the  matched  transition  of  the  coaxial  exciter  was  constructed. 
Of  course,  the  efficiency  of  the  source  is  independent  of  the  feeding 
structure,  but  substantial  coupling  of  power  from  the  closed  waveguide 
to  the  surface  waveguide  is  obtained  only  when  a  good  impedance  match 
is  provided  at  the  aperture  plane.   Otherwise,  a  very  high  standing 
wave  ratio  exists  in  the  vicinity  of  the  aperture  and  most  of  the 
power  delivered  by  the  feeding  waveguide  is  dissipated  as  loss  at  this 
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position   A  matched  transition  is  of  utmost  importance  when  one  wishes 
to  excite  any  type  of  surface  waveguide. 

Finally,  good  experimental  verification  of  the  pole,  saddle  point 
method  of  calculating  efficiency  was  obtained., 
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APPENDIX  A 

In  order  to  solve  for  the  constant  Ds  we  write  the  system  of 
equations  (3,12)  in  matrix  form: 


\l 

*L2 

^3 

0 

A 

*» 

\2 

^3 

0 

B 

0 

^2 

^3 

*34 

C 

0 

A42 

A43 

A44 

D 

ioe   b 
o 

0 
0 


(A.l) 


where 


\l   =   "  (Vlb)  Jl   (Via) 

*L2  =  "  *ll 

^3  =  (Vlb)  Yl  (Via) 


*2l  =  VVia) 
^2  =  "Al 
^23="  Yl(Via) 


^2  = 

^(Vjb) 

^3  = 

Yl(v1b) 

^4   = 

(v  b) 
o 

A42  = 

(vlb)    j; 

O^b) 

A43  =  (vlb)  Yx  (vlb) 

A44  "  VV»  Hia)(Vob)+(er-l)H1(1>(Vob) 


Denote  the  coefficient  matrix  of  (Asl)  by  j~A.  .  "1  and  let  the  cof actor 
A.  .  be  designated  X.  ..   The  solution  of  (A.l)  is 


of 


Det 


[AiJ 


xll 

X2! 

X31 

X41 

X12 

X22 

X32 

X42 

*13 

X23 

X33 

X43 

\4 

X24 

X34 

X44 

-  ia)€  b 
0 
0 
0 


(A,  2) 


which  reduces   to 


HI 


•Mb 


Det 


[Au] 


r~                 -1 

\l 

X12 

\3 

\4 

Evaluating  Det    ["A.  .1    ,    we  obtain 

DGt    [Aij]       =  *L1  (A22  *23  A44  +  ^23  A»4  A42> 

-  *11  (A32  ^3  A44  +  ^2  ^4  A43) 

-\2  (A21  ^3  A44  "  ^1  ^4  A43) 

+  ^3  (A21  ^2  A44  "  *2l  ^4  A42) 


(A3) 


(A, 4) 


since  A   =  -  AL2  and  Ag  =  -  \   ,       (A.  4)  reduces  to 


DGt  [Aij]  "  (A13  ^1  "  *ll  A23)  (A32  A44  "  ^4  A42>' 


(A.  5) 


We  note  that 


(AL3  *21  "  *L1  A23)  =  (Vlb)  [Jl(Via)  Yl  (Via)  "  Yl(Via)  Jl  (Vla)J 
=  (Vlb)^aT 


(A.  6) 


since  the  bracketed  term  is  the  Wronskian  of  J  (v  a) ,  Y  (v  a)  which  is 
|  equal  to  [2/7T  (v  a)].   Substituting  (A ,6)  into  (A,  5)  we  obtain 


2b 


Det  [Aij]    ■  15    (A32  A44   "  *34  A42> 


(A.7) 


The  cof actor  X        is  given  by 


From   (A, 3),     (A. 7)    and    (A. 8)    we   have 


iJJ  ff    (A32  A44   "  ^4  A42> 


Consider 
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14  =  "  Aai  (A32A43  -^a^a>-  <A-8> 


•    *   h  ^  *  A21(A32A43   -  A42A33) 

D  =  "    1Cj€lb  Bet    FA.  .1    =  +    laCl*      2    ,A        A T" — — ■  <A'9> 

L   ijJ  =   (A  „  A_   -  A_  A     ) 


(A32  A43   "  A42  ^3)    =    (Vlb)    LVV*    Tl    (V>    "   \  <V>    V   <V>] 

=  (vib)^bT  =  l  (A'10) 

because  of  Wronskian  relation.      Utilizing  Bessel   function  recurrence 
relations  we  determine   that 


(A32  A44   "  ^4  A42>    =    (Vlb)    Jo(Vlb)    *i™    (Vob)    "  V\>b)  Jl(vlb)Ho(1)  <Vob) 


(A.1D 


Substituting    (A. 10)    and    (A, 11)    into    (A,  9)    yields 


+   ±cj€   a  J   (v  a) 

>- — i_L_i (A. 12) 

(vib)   Jo(Vlb)   Hx(   >(v  b)    -  €r(vob)    Jl(Vlb)   H/        (vob) 
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APPENDIX  B 

The  graphical  solution  of  Eqs .  (4„9)  and  (4.10)  is  illustrated  in 
Fig.  4  of  Chapter  4.   Equations  (4.9)  and  (4.10)  are 


JQ(X)         KQ(£) 

T[ao  '  £r(e>  K]T¥) 


=   0 


(4.9) 


r2  =  (knb)2  (s  -  i)  =  r2 

u  r 


(4  JO) 


where  we  omit  the  subscript  1  of  X  for  convenience  of  notation. 

Referring  to  Fig,  4,  we  see  that  the  intersection  of  the  circle  and 
the  first  branch  of  X  =  f(£)  is  the  unique  solution  of  (4.9)  and  (4.10)0 
Consider  the  magnified  view  at  the  point  of  intersection  of  the  curves, 
as  shown  in  Fig.  21. 


FIGURE  21.   SOLUTION  OF  THE  MODE  EQUATION 
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Denote  points  on  the  circle  corresponding  to  values  of  %   by  X,  where 
X  =  sj  R2  -  i2   from  (4.10).   The  points  X,  which  define  the  curve  with 
positive  slope,  are  solutions  of  (4.9),  which  may  be  written  in  the  form 


J0(X) 
-Xj7W-V  =  °  (B.l) 

k0(6) 

where  v  =  £   (£)  _.  /gv   is  a  constant  for  any  selected  £  and  value  of  %. 
»     r     K  (5)  r 

We  see  that  the  points  X  are  simply  real-roots  of  the  equation 

J0(X) 
-  X   (  v  -  y  =  0.   Let  the  point  of  intersection,  which  is  the  exact 

solution  of  (4.9)  and  (4.10),  have  the  coordinates  %     and  X  .   The  basic 

s      s 

method  for  determining  the  point  of  intersection  using  the  digital 

computer  is  illustrated  in  Fig.  21.   It  is  evident  that  for  any  value 

£  less  than  §  ,  X'  >  X,  and  when  I  is  greater  than  I  ,  x'<X.   By 

computing  X'  and  X  for  sufficiently  small  increments  of  £,  it  is  possible 

to  determine  the  neighborhood  of  intersection  very  accurately. 

Two  programs  were  written  to  determine  t   and  X  to  an  accuracy  of 

s      s 

18 
four  decimal  places.   The  British  Association  Mathematical  Tables 

tabulate  the  functions  K  (I)  and  K  (6)  for  .01  increments  of  I.   The  first 

program  used  the  tabulated  functions  to  determine  the  intersection  as 

occurring  between  §.  and  I.  ,,  where  the  interval  %.    _  -  I.  =0.01.   The 
i      l+l  l+l    i 

second  program  included  an  interpolation  routine  which  subdivided  the 
previously  determined  intersection  interval  into  100  parts  and  thereby 

yielded  §  and  X  accurate  to  four  decimal  places.   A  description  of 

s      s 

each  program  is  given  in  Sections  B.l  and  B.2,  respectively. 
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B.l 

The  ILLIAC  registers  are  limited  to  numbers  in  the  range  -1  <  p  <  1. 
All  numerical  computations  must  yield  results  which  are  less  than  one  in 
absolute  value;  otherwise  the  registers  overflow  and  one  obtains  meaning- 
less results.   Fractions  read  into  the  computer  must  also  be  in  the  range 
-1  <  p  <  1 .   Scaling  of  computations  is  almost  always  necessary  and  this 
program  was  no  exception.   The  program  was  written  using  the  scaled 

quantities  k^b/8,  £  /4,  £/4   and  X/4.   If  the  functions  K„(£)  and  KAt) 
u      r  U         1 

exceeded  unity,  the  values  actually  read  into  the  computer  were  1/2  K_(£) 

and  1/2  K, (I) .   Selecting  8-  =  2.56  and  a  particular  value  of  k_b/8 ,  the 
1  r  L) 

routine  computed  ten  values  of  £/4,  X/4,  and  X'/4  for  .01  increments  of  5 
in  the  vicinity  of  the  point  of  intersection.   The  program  utilized  the 
following  ILLIAC  library  routines: 

HI  -  Inverse  Interpolation;  finds  a  real  root  of  f(X)  =  0  between 
two  values  X  and  X  of  the  argument  of  f (X) „ 
N12  -  Infraput;  a  routine  to  read  in  up  to  12  digit  fractions  or 

integers  which  are  punched  on  tape.   The  sequence  of  fractions 
or  integers  is  placed  in  memory  locations  n,  n+1,  n+2,  .... 
P16  -  Infraprint;  this  integer,  fraction  print  routine  prints  the 
contents  of  the  accumulator  up  to  12  decimal  places,  with 
optional  choice  of  decimal  point,  sign,  etc. 
Rl  -  Square  Root;  uses  Newton's  method  to  replace  p  in  the  accumur 

lator  by   I p, 
V8  -  Ordinary  Bessel  Functions;  computes  J  (X)  ,  J^X),  Y  (X)  ,  and 
Y  (X)  to  an  error  of  less  than  2    ,  using  series  expansions 
for  the  functions. 
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Library  routine  HI  requires  the  user  to  write  an  auxiliary  subroutine 

which  computes  f (X) .   The  subroutine  which  was  written  refers  to  V8  in 

JQ(X) 

order  to  compute  f(X)  =  -  X    .    .  -  v, 

J  V.X.J 

The  program  commenced  by  reading  in  ten  values  of  K_(6)  and  placing 
them  in  memory  locations  500  through  509;  ten  corresponding  values  of 
K  (§)  were  placed  in  locations  510  through  519,  while  placing  Initial 
6/4  in  520,  X^/4  in  521,  X2/4  in  522,  and  k  b/8  in  523.   Note  that  X 
and  X  are  arbitrary  bounds  of  X  within  which  a  root  of  f(X)  =  0  is 
assumed  to  occur.   As  an  example,  consider  the  case  k  b  =  3.0,  which 
yields  R  =  3.75.   Figure  4  indicates  the  intersection  at  approximately 
6  =  1.97,  X  =  3.2,  so  6  was  chosen  to  have  the  values  6  =  1,92,  193, 
1.94,  ....,  2.01  with  X  =  3.16  and  X  =  3.26.   Starting  with  the  initial 

£r  6  Ko(l) 

value  of  5/4,  1/1 6v  =  -j—  T  v    }£\      was  computed,  and  library  routine  HI 
was  entered  to  obtain  a  scaled  root  X/4  of  the  scaled  equation 

-  Yg  (X)  JQ(X)/J1(X)  -  ~Y  =  0.   Using  kQb/8,  £r/4,  and  1/4,  X'/4  was 


calculated  from  the  equation  X'/4  =  -jJCk  b)   (£  -1)  -  6  .   The  quantities 

k  b/8,  6/4,  X/4,  and  X9/4  were  printed  to  four  decimal  places,  6  was 

increased  by  .01  and  the  computation  was  repeated.   The  cycle  continued 

until  all  ten  values  of  6  had  been  used,  which  completed  the  computations 

for  that  particular  k  b„   New  constants  for  a  different  k  b  were  read  from 

the  tape  and  the  entire  routine  was  repeated.   A  first  approximation  of  6 

and  X  for  the  six  values  of  k^b  was  obtained  in  this  manner.   The  results 
s  0 

are  presented  in  Table  B.l.   Preceding  each  set  of  data  is  the  value  of 
k Qb/8.   The  first,  second,  and  third  columns  give  6/4,  X/4,  and  X'/4, 
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respectively.   The  arrow  indicates  the  interval  of  £  within  which  inter- 
section occurs.   The  tape  format  and  order  codes  of  the  routines  are 
as  follows: 
Tape  Format 

00  20  K 
(N12)  00  K 
(P16)  00     K 

00   3  K 

00   F  00   600  F 

00  115  K 
(Rl)   00   K 

00   K 

00  600  K 

00  525  K 

24  525  N 


50 

500  F 

50 

L 

26 

20  F 

L5 

520  F 

L4 

(8) 

40 

(9) 

92 

707  F 

92 

131  F 

92 

515  F 

Routine  N  12 
Routine  P  16 

pre-set  parameter  for  HI 

Routine  HI 

Routine  Rl 

Routine  V8 

Auxiliary  Subroutine 

Main  Routine 

Stop  -  Transfer  control  to  main  routine 

Main  Routine 
read  in  24  fractions 


.088 


form  end  test  constant 


L5 

523  F 

5 

J4 

3  F 

50 

5  L 

6 

26 

59  F 

L5 

520  F 

7 

40 

(1) 

)92 

131  F 

8 

92 

515  F 

L5 

(1) 

9 

J4 

4  F 

50 

9  L 

10 

26 

59  F 

print  kQb/8 


initial   6/4  in   (1) 


92 


975   F 


print  £./4 


4  spaces 


11* 

50 

(1) 

7J 

500  F 

12* 

50 

(2) 

66 

510  F 

13 

7J 

(3) 

40 

(4) 

14 

L5 

521  F 

40 

10  F 

15 

L5 

522  F 

40 

12  F 

1  Jl  , 


4  rinQ 


16  V  in   (4> 


4  Xx   in  10 


4  X2  in  12 


16 


17 


18 


19 


20 


21 


22 


23 


24 


25 


26 


27 


28 


JO 

(2) 

50 

16L 

26 

115  F 

JO 

(2) 

J4 

4  F 

50 

18  L 

26 

59  F 

92 


975  F 


L5 

(3) 

LO 

(5) 

40 

(6) 

50 

523  F 

7J 

523  F 

40 

(7) 

50 

(6) 

7J 

(7) 

40 

(7) 

L5 

(1) 

10 

2  F 

40 

(6) 

50 

(6) 

7J 

(6) 

40 

(6) 

L5 

(7) 

LO 

(6) 

waste 
enter  HI 


waste 


print  X^/4 


4  spaces 


4  <£,-!)    in  (6) 


(kb/8)2  in  (7) 


(^[(kb)^  (£-!)]  m  (7) 


16   in  (6) 


1  n2 


(— >*  [(kb)2  (^-1)  -  e.2]  in  A 


90 


50 

28  L 

29 

26 

149  F 

50 

(2) 

30 

00 

2  F 

JO 

(2) 

31 

J4 

4  F 

50 

31  L 

32 

26 

59  F 

L5 

(1) 

33 

LO 

(9) 

<= 

36 

39  L 

34 

L5 

(1) 

L4 

(10) 

35 

40 

(1) 

50 

11  L 

36 

K5 

F 

42 

11  L 

37 

50 

12  L 

K5 

F 

38 

42 

12  L 

enter  Rl ,  leave  r-r  X .  '  in  A 
16   l 


clear  Q 


form  i  X,  ' 

4   i 


waste 


print  —  X. ' 


i.         |  +  .088 
test:   (^  -  -^ ) 


advance  £. 


i.    +    .01 

l 


in  (1) 


advance  (500)  address 


advance  (510)  address 


22 


7  L 


transfer  control  to  7  L 


39 


L5 
42 


(11) 
11  L 


reset  (500) 
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40      L5 

(12) 

42 

12  L 

reset  (510) 

41      59 
LO 

(2) 
523  F 

—  in  A,  clear  Q 

test:  .    _ 

4.0  -  kb 

8 

42  < 36 

L 

OF 

F 

43  (1)  00 

F 

00 

F 

44  (2)  00 

F 

00 

F 

45  (3)  40 

F 

00 

140  000  000  000  J 

46  (4)  00 

F 

00 

F 

47  (5)  00 

F 

00 

250  000  000  000  J 

48  (6)  00 

F 

00 

F 

49  (7)  00 

F 

00 

F 

50  (3)  00 

F 

00 

22  000  000  000  J 

51  (9)  00 

F 

00 

F 

52  GO)  00 

F 

00 

2  500 

000  000  J 

92 


(11) 

00 

F 

00 

500 

F 

(12) 

00 

F 

00 

510 

F 

Auxiliary  Subroutine 

40 
K5 

18 

F 
F 

rescue  X. /4 

42 

8 

L 

plant  link 

L5 

18 

F 

50 

2 

F 

enter  V8 ,  leave  JQ(X)  •  5 

50 

2 

L 

-19 
and  J  (X)  -2    in7. 

26 

158 

F 

L5 

18 

F 

10 

2 

F 

form  X/16,  store  in  19. 

40 

19 

F 

50 

6 

F 

79 

19 

F 

50 

(2) 

66 

7 

F 

1    Jo(x)   1 

f°rm  '  16  X  J]L(X)  "  16  T 

S5 

F 

L0 

(4) 

JO 

(2) 

waste 

22 

F 

•19 
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TABLE  B.l 

275 

1375 

.6711 

.6730 

1400 

.6725 

.6725  -<s 

1425 

.6739 

.6720 

1450 

.6753 

.6715 

1475 

.6767 

.6709 

1500 

.6781 

.6704 

1525 

.6795 

.6698 

1550 

.6809 

.6692 

1575 

.6823 

.6687 

1600 

.6837 

.6681 

325 

3075 

.7508 

.7514 
.7503 

3100 

.7517 

3125 

.7526 

.7493 

3150 

.7535 

.7482 

3175 

.7544 

.7472 

3200 

.7553 

.7461 

3225 

.7562 

.7450 

3250 

.7571 

.7440 

3275 

.7580 

.7429 

3300 

.7588 

.7418 

375 

4800 

.8013 

.8044 

4825 

.8018 

.8029 

4850 

.8024 

.8014  "^ 

4875 

.8030 

.7999 

4900 

.8035 

.7984 

4925 

.8041 

.7968 

4950 

.8046 

.7953 

4875 

.8052 

.7937 

5000 

.8057 

.7921 

5025 

.8063 

.7906 

425 

6525 

.8335 

.8375 

6550 

.8339 

.8355 
.8335 

6575 

.8343 

6600 

.8346 

.8316 

6625 

.8350 

.8296 

6650 

.8354 

.8276 

6675 

.8357 

.8256 

6700 

.8361 

.8235 

6725 

.8365 

.8215 

6750 

.8368 

.8194 
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TABLE  B.l  (CONTINUED) 


8200 

.8548 

.8576 

8225 

.8551 

.8552 
.8528 

8250 

.8554 

8275 

.8556 

.8504 

8300 

.8559 

.8479 

8325 

.8561 

.8455 

8350 

.8564 

.8430 

8375 

.8567 

.8405 

8400 

.8569 

.8380 

8425 

.8572 

.8355 

525 

9750 

.8692 

.8771 

9775 

.8694 

.8743 

9800 

.8696 

.8715 
.8687 

9825 

.8698 

9850 

.8700 

.8658 

9875 

.8702 

.8630 

9900 

.8704 

.8601 

9925 

.8705 

.8572 

9950 

.8707 

.8543 

9975 

.8709 

.8514 

A  parameter  tape  was  placed  in  the  photoelectric  tape  reader  after 

the  main  program  tape  had  been  read  into  the  computer.   The  parameter 

tape  listed  the  Bessel  functions  and  constants  which  were  required  for 

the  computations.   The  fractions  had  to  appear  on  the  tape  in  the  proper 

sequence,  which  was:   Ten  values  of  Kft  (!)  ,  ten  corresponding  values  of 

K  (!)  ,  and  initial  J?/4,  X  /4,  X  /4,  and  k  b/8.   The  total  parameter  tape 

was  comprised  of  six  sets  of  24  fractions  each,  where  24  fractions  were 

required  for  each  value  of  k  b.   The  entries  on  the  tape  for  k  b  =  3.0 

o  o 

were  as  follows: 


+12569493 
+12415057 
+12262656 
+12112263 
+11963846 
+11817379 
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+11672831 

+11530177 

+11389387 

+112504361 

+15546507 

+15341358 

+15139211 

+14940014 

+14743718 

+14550275 

+14359637 

+14171756 

+13986588 

+138040877 

+48 

+  79 

+  815 

+  375  N 


A  decimal  point  before  the  first  digit  is  assumed  by  the  ILLIAC.   The 
first  ten  entries  are  K  (1.92),  K  (1,93),  .  ,..,  K  (2.01)   and  the  second 
ten  fractions  are  K  ( 1.92) ,  K  (1.93),  ....,  K  (2.01).   The  terminating 
symbol  N  following  k  b/8  tells  ILLIAC  to  stop  reading  in  fractions  and 
control  is  returned  to  the  main  routine. 
B.2 

The  second  program  utilized  library  routine  II  in  addition  to  all 
of  the  library  routines  which  were  used  in  the  first  program.   Routine 
II  computes  an  interpolate  from  a  table  of  function  values  stored  in  the 
computer  memory.   It  uses  Neville's  method  of  successive  linear  inter- 
polations.  The  user  specifies  the  depth,  d,  of  interpolation,  which 

1  implies  the  highest  order  "difference"  used  in  the  interpolation  formula. 

1  o  t      K0(^ 
i  In  this  case,  the  depth  was  specified  as  four.   The  function  —  5  (5)    (£\ 

is  fairly  linear  and  more  amenable  to  interpolation  than  the  K  Bessel 

function.   Consequently,  values  of  this  function  were  read  into  the  computer 

for  t   in  the  neighborhood  of  £  and  X  ,  and  interpolation  was  performed 
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within  the  sequential  table  of  values.  The  first  program  had  determined 
the  intersection  as  occurring  between  5  and  §.  ,,  where  5...  -  £,  -    .01 


i+1 


i+1 


The  second  program  increased  s  within  this  interval  in  increments  of 

i         ve> 

.0001,  interpolating  to  obtain  —  £   (£)    (r      for  each  £,  until  inter- 
section was  detected  by  the  method  shown  in  Fig.  21.   When  intersection 
was  detected,  the  computer  printed  out  %,   X,  and  X'  on  each  side  of 
crossover  in  the  following  form: 

k  b/8 
o 


K 


hi 


4  i+1 


K 


I  x 

4  *i+l 


ix'. 

4   l 


—    Y> 

4  A  i+1 


Table  B.2  presents  the  numerical  results  for  each  of  the  six  values  of 

k  b.   Table  I  was  constructed  from  these  data, 
o 


TABLE  B.2 


14005 

.67251 

.67252 

14007 

.67252 

.67252 

325 

30820 

.75107 

.75107 

30822 

.75108 

.75106 

375 

48380 

.80214 

.80215 

48382 

.80214 

.80213 

425 

65672 

.83415 

.83415 

65675 

.83415 

.83413 

475 

82260 

.85510 

.85512 

82262 

.85510 

.85510 

525 

98157 

.86969 

.86972 

98160 

.86969 

.86969 

97 


The  procedure  for  calculating  X  and  X"  was  similar  to  that  used  in 
the  first  program.   The  only  essential  difference  was  the  addition  of 
routine  II  to  perform  the  interpolation.   The  tape  format  and  order  codes 
written  for  the  main  routine  and  auxiliary  subroutine  are  as  follows: 
Tape  Format 


Routine  N12 
Routine  PI 6 
Routine  Rl 


prerset  parameter  for  HI 
Routine  HI 

pre-set  parameter  for  II 

Routine  II 

Routine  V8 

Main  Routine 

Auxiliary  Subroutine1 

Stop  -  Transfer  Control  to  Main  Routine 
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Main  Routine 


50 

21  F 

50 

L 

26 

200  F 

JO 

(1) 

50 

54  F 

50 

2  L 

26 

200  F 

26 

75  L 

-3-50     190  F 


50 


4  L 


26 

200  F 

JO 

(1) 

50 

40  F 

50 

6  L 

26 

200  F 

JO 

(1) 

52 

45  F 

50 

8  L 

26 

200  F 

read  in  12  fractions,  locations  21 
through  32 

waste 

read  12  fractions  into  54  through  65 


read  6  fractions  into  190  through  195 


waste 


read  parameter  constants  into  40  through  44 


waste 


-39 
read  integer  a  *  2    into  45 


gg 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


21 


L5 

42 

92 

92 

92 

JO 

L5 

J4 

50 

26 

L5 

40 

L5 

50 

50 

26 

00 

00 

00 

40 

L5 

40       10  F 

L5       42  F 

40       12  F 


45  F 

17  L 

707  F 

131  F 

515  F 

(1) 
43  F 

3  F 
13  L 

239  F 

40  F 
(2) 
(2) 

4  F 
16  L 

338  F 

F 

F 
10  000  000  000  J 
(3) 

41  F 


insert  a  into  II  entry 

number  shift 

carriage  return  and  line  feed 

delay 

waste 

enter  PI 6,  print  kb/8 


puts  initial  £   into  (2) 


current  £   in  A 


enter  II,  interpolates  to  leave 


TiY  inA 


store  interpolate  in  (3) 


X  /4  in  10 


X2/4  in  12 


100 


22 


23 


24 


25 


26 


27 


28 


29 


30 


31 


32 


33 


34 


50 

21  L 

26 

304  F 

40 

(4) 

L5 

(2) 

10 

2  F 

L4 

44  F 

40 

(5) 

50 

43  F 

7J 

43  F 

40 

(6) 

50 

(6) 

7J 

(7) 

40 

(6) 

L5 

(5) 

10 

2  F 

40 

(8) 

50 

(8) 

7J 

(8) 

40 

(8) 

41 

F 

L5 

(6) 

L0 

(8) 

50 

32  L 

26 

295  F 

50 

(1) 

00 

2  F 

40 

(9) 

enter  HI,    leave   root   X. /4   in  A 

l 

store  X. /4   in   (4) 


forms   £ .  /4   from  t 
l  t 


store   £;/4   in   (5) 


(kb/8)      in   (6) 


(— )2     [(kb)2(£fl)]    in   (6) 


§./16   in   (8) 


(i/16)2   in   (8) 


clear  A  and  location  0 


(T6")2   t<kb>2  <£f  D  -  ^2]  in  A 


enter  Rl ,  leave  —  X. J  in  A 
16   i 

clear  Q 


X^/4  in  (9) 


101 


35 


36 


37 


38 


39 


40 


41 


42 


43 


44 


45 


46 


47 


48 


L5 

(4) 

L0 

(9) 

-32 

41  L 

L5 

(5) 

40 

(10) 

L5 

(4) 

40 

(ID 

L5 

(9) 

40 

(12) 

L5 

(2) 

L4 

(13) 

40 

(2) 

22 

15  L 

92 

131  F 

92 

515  F 

L5 

(10) 

J4 

5  F 

50 

43  L 

26 

239  F 

92 

975  F 

JO 

(1) 

L5 

(ID 

J4 

5  F 

50 

46  L 

26 

239  F 

92 

975  F 

JO 

(1) 

test:   X./4  -  X.  /4 
i      i 

transfer  to  print 


§  /4  in  (10) 


X./4  in  (11) 


X.  /4  in  (12) 


advance  6  by  .0001 


carriage  return,  line  feed,  delay 


print  £ . /4 


four  spaces 
waste 


print  X./4 


four  spaces 
waste 
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49 


50 


51 


52 


53 


54 


55 


56 


57 


58 


60 


61 


L5 

(12) 

J4 

5  F 

50 

49  L 

26 

239  F 

92 

131  F 

92 

515  F 

L5 

(5) 

J4 

5  F 

50 

52  L 

26 

239  F 

92 

975  F 

JO 

(1) 

L5 

(4) 

J4 

5  F 

50 

55  L 

26 

239  F 

92 

975  F 

JO 

(1) 

L5 

(9) 

J4 

5  F 

50 

58  L 

26 

239  F 

59 

(1) 

LO 

43  F 

36 

6  L 

OF 

F 

00 

F 

print  X  /4 


carriage  return  and  line  feed 
delay 


print  i±+1/4 


waste 


print  X.  ,  /4 
l+l 


print  X.  +  1/4 


in  A,  clear  Q 


test: 


4.0  -  kb 
8 


transfer  control  to  read  in  new  parameters 
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62  (1)   00  F  00  F 

63  (2)   00  F  00  F 

64  (3)   00  F  00  F 

65  (4)   00  F  00  F 

66  (5)   00  F  00  F 

67  (6)   00  F  00  F 

68  (7)   00  F  00   390  000  000  000  J 

69  (8)   00  F  00  F 

70  (9)   00  F  00  F 

71  (10)  00  F  00  F 

72  (11)  00  F  00  F 

73  (12)  00  F  00  F 

74  (13)  00  F  00       100  000  000  J 


75  50  91  F 
50   75  L 

76  26  200  F 

-*—  26   4  L 


read  6  fractions  into  91  through  96 


40 

18  F 

K5 

F 

42 

8  L 

L5 

18  F 

50 

2  F 

50 

2  L 

26 

389  F 

Auxilliary  Subroutine 
rescue  X. /4 

plant  link 
X/4  in  A 


enter  V8< 


J0(X) 


■JX(X) 


2~19  in  6 


2-19  in  7 
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L5 

18  F 

10 

2  F 

40 

19  F 

50 

6  F 

79 

19  F 

50 

(1) 

66 

7  F 

S5 

F 

L0 

(3) 

JO 

(1) 

-< 22 

F 

form  XA6,  store  in  19  F 


TS  Vx>  '  2"19  in  A 


clear  Q 

X  Jo(x) 

16  ^(Xj.) 


in  A 


X  J0<X)    1  v   .   A 
16VX)-IiY  lnA 


waste 


The  parameter  tape  which  was  read  into  the  computer  after  the  main 

i        ve> 

program  tape  listed  the  function  values,  —  £   (I)  T-  , g x  ,  required  by 

lb  r     K  \b) 

the  interpolation  routine.   It  also  included  constants  telling  the 
computer  where  the  table  was  stored,  the  initial  value  of  b,  X  /4,  X  /4, 

1      2i 

and  k  b/8.   The  parameter  tape  as  it  was  read  into  ILLIAC  is  as  follows: 
o 

+14217  205 
+14363  079 
+14509  095 
+14655  256 
+14801  554 
+14947  990 
+45745  8333 
+45902  4109 
+46059  0048 
+46215  6141 
+46372  2411 
+46528  8834N 


+05055  518 
+05089  839 
+05216  642 
+05343  995 
+05471  883 
+05600  291 
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+35300  6528 

+35455  7627 

+35610  9008 

+35766  0675 

+35921  2618 

+36076  4848N 

+24685  093 
+24837  399 
+24989  766 
+25142  182 
+25294  660 
+25447  191N 

+55638  3502 
+55795  7729 
+55953  2083 
+56110  6549 
+56268  1103 
+56425  5788N 

+  555+  66+  68+275+  HfN 

+  23+74+76+ 325+25N+N 

+93+79+81+ 375+ 25N+N 

+  62+82+84+  425+5N+N 

+  285+84+86+475+75N+N 

+92+86+88+525+75N+100N 
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APPENDIX  C 

The   integral    (6.6)   which  was  needed   in  order   to  compute   the 
excitation  efficiency   is 


7T/2 


cos   9  F      (9)  d9  (6.6) 


where 


F(9)   =  


J,  (k  a  w) 
1      o 


w  J   (k  bw)   H. (1)(k  b  cos  9)   -  €      cos  9  J   (k  bw)   H   (1)(k  b   cos  9) 
ool  o  r  looo 


and 


w  = 


€      -   sin2 
r 


The  program  which  was  written  to  determine  N  utilized  ILLIAC 

library  routine  E2,  which  computes  an  approximation  to  the  integra] 

b 
7—     f  (x)  dx  by  Simpson'  s  rule.   This  routine  requires  that  the 


user  supply  tabulated  values  of  the  integrand  at  an  even  number  of 

equally  spaced  intervals.    The  range  of  integration,  which  was  zero 

to  88  degrees,  was  subdivided  into  44  intervals  of  two  degrees  each. 

For  any  k  a  and  k  b,  the  integrand  was  computed  at  the  45  points 
o       o 

0,  2°,  4°,  6°,  ...  ,  86°,  88°.   Since  the  integrand  contains  F(9) , 
which  involves  a  division,  it  was  necessary  to  scale  the  numerator 
whenever  it  exceeded  the  denominator  in  order  to  maintain  |  F(9)  |  <  1; 
otherwise  the  capacity  of  the  ILLIAC  registers  would  be  exceeded.   This 
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was  done  by  multiplying  the  numerator  by  1/2  as  often  as  needed  until 
all  45  values  of  F(0)  were  within  range.   Consequently,  the  computed 
integral  was  scaled  by  an  integer  number  of  multiplications  by  the 
factor  1/2.   The  scaled  integral,  I,  which  was  computed  was  equal  to 

88  £_ 
n90  2      p+l 

I  =  |8^     J  <2>      C°S  9   I  F(9)  I 

90  2 

where  p  is  the  integer  scale  factor. 

R 
It  follows  from  (6.6)  and  (Col)  that  N  is  obtained  from  I  by 


For  a  particular  choice  of  k  b,  the  integral  (CI)  was  calculated 

o 

for  k  a  equal  to  0.2,  0.4,  0.6,  ...  until  k  a  took  on  the  value  k  b; 
o  o  o 

then  k  b  was  changed  to  a  new  value,  and  the  process  was  repeated.   The 
o 

scaled  quantities  k  b/8  and  k  a/8  were  used  in  the  computations.   The 
o         o 

results  were  printed  in  four  columns  which  listed,  respectively, 

k  b/8,  k  a/8,  p,  and  I.   Relation   (C.2)  was  then  used  to  compute  the 
o      o 


(C2) 


NR 


corresponding  value  of  N  .   Table  CI  presents  the  results  as  printed  out 

by  ILLIAC  for  k  b/8  equal  to  0.425,  that  is,  k  b  equal  to  3.4. 
o  o 

The  program  which  was  written  to  compute  the  integral  (C.l)  was 
comprised  of  143  order  pairs.   Due  to  its  length,  the  order  code  will 
not  be  given  in  detail.   However,  it  should  be  mentioned  that  library 
routines  V8  (Bessel  functions),  T5  (sine-cosine),  Rl  (square  root), 
P16  (print  routine),  and  E2  (integration)  were  utilized  by  the  main 
program. 


108 


TABLE  CI. 

k  b/8 
o 

k  a/8 
o 

P 

I 

.425 

.025 

1 

.0304 

.425 

.050 

1 

.1136 

.425 

.075 

1 

.2280 

.425 

.100 

2 

.1722 

.425 

.125 

2 

.2167 

.425 

.150 

3 

.1184 

.425 

.175 

3 

.1142 

.425 

.200 

3 

.0974 

.425 

.225 

2 

.1450 

.425 

.250 

2 

.0913 

.425 

.275 

1 

.0914 

.425 

.300 

1 

.0334 

.425 

.325 

1 

.0149 

.425 

.350 

1 

.0310 

.425 

.375 

1 

.0682 

.425 

.400 

1 

,1097 

.425 

.425 

1 

.1403 
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ATTN:  Dr,  A.  Alford 
299  Atlantic  Ave. 
Boston  10,  Massachusetts 
M/F  Contract  AF33(038)-23700 

Bell  Aircraft  Corporation 
Attn:   Mr.  J.D,  Shantz 
Buffalo  5,  New  York 
M/F  Contract  W-33(038)-l4l69 

Chief 

Bureau  of  Ships 

Department  of  the  Navy 

ATTN:   Code  838D,  L.E.  Shoemaker 

Washington  25,  D.C. 

Director 

Naval  Research  Laboratory 

ATTN:   Dr.  J.I.  Bohnert 

Anacostia 

Washington  25,  D.C. 
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National    Bureau   of   Standards 
Department    of   Commerce 
ATTN:       Dr.    A,G.    McHish 
Washington  25,    D.C. 

Director 

U.S.   Navy  Electronics   Lab. 

Point   Loma 

San  Diego   52,    California 

Commander 

USA  White  Sands   Signal   Agency 

White  Sands   Proving  Command 

ATTN:      SIGWS-FC-02 

White  Sands,    N.M. 

Consolidated  Vultee  Aircraft  Corp, 

Fort  Worth  Division 

ATTN:      CR.    Curnutt 

Fort  Worth,    Texas 

M/F  Contract  AF33(038)-21117 

Dome   &  Margolin 

29  New  York  Ave. 

Westbury 

Long   Islai.d,    New  York 

M/F  Contract  AF33( 616) -2037 

Douglas  Aircraft   Company,    Inc. 

Long  Beach  Plant 

ATTN:      J.C.    Buckwalter 

Long  Beach  1,    California 

M/F  Contract  AF33( 600) -25669 

Boeing  Airplane  Company 

VTTN:      F.    Bushman 

7755  Marginal   Way 

Seattle,    Washington 

i/F  Contract  AF33(038)-31096 

'hance-Vought  Aircraft   Division 
Inited  Aircraft   Corporation 
iTTN:      Mr.    F.N.    Kickerman 
'HRU:      BuAer  Representative 
•alias,    Texas 

onsolidated-Vultee  Aircraft  Corp, 
TTN:      Mr.    R.E.    Honer 
.0.   Box  1950 

ian  Diego  12,   California 
/F  Contract  AF33( 600) -26530 


Grumman  Aircraft   Engineering  Corp. 
ATTN:       J.S.    Erickson, 

Asst .    Chief,    Avionics    Dept . 
Bethpage 

Long    Island,    New  York 
M/F  Contract   NOa(s)    51-118 

Hallicrafters   Corporation 

ATTN:      Norman   Foot 

440  W.    5th  Avenue 

Chicago,    Illinois 

M/F  Contract   AF33(600)-26117 

Hoffman  Laboratories,  Inc. 

ATTN:   S.  Varian 

Los  Angeles,  California 

M/F  Contract  AF33<600)-17529 

Hughes  Aircraft  Corporation 

Division  of  Hughes  Tool  Company 

ATTN:   D.  Adcock 

Florence  Avenue  at  Teale 

Culver  City,  California 

M/F  Contract  AF33(600)-27615 

Illinois,  University  of 

Head,  Department  of  Electrical  Engineering 

ATTN:   Dr.  E.C.  Jordan 

Urbana,  Illinois 

Johns  Hopkins  University 
Radiation  Laboratory 
ATTN:   Librarian 
1315  St.  Paul  Street 
Baltimore  2,  Maryland 
M/F  Contract  AF(33)616)-68 

Glenn  L.  Martin  Company 

Baltimore  3,  Maryland 

M/F  Contract  AF33( 600) -21703 

McDonnell  Aircraft  Corporation 
ATTN:   Engineering  Library 
Lambert  Municipal  Airport 
St.  Louis  21,  Missouri 
M/F  Contract  AF33(600)-3743 

Michigan,  University  of 
Aeronautical  Research  Center 
ATTN:   Dr.  L.  Cutrona 
Willow  Run  Airport 
Ypsilanti ,  Michigan 
M/F  Contract  AF33(038)-21573 
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Massachusetts  Institute  of  Tech. 
ATTN:   Prof.  H.J.  Zimmermann 

Research  Lab.  of  Electronics 
Cambridge,  Massachusetts 
M/F  Contract  AF33(616)-2107 

North  American  Aviation,  Inc. 
Aerophysics  Laboratory 
ATTN:   Dr.  J. A.  Marsh 
12214  Lakewood  Boulevard 
Downey,  California 
M/F  Contract  AF33( 038) 18319 

North  American  Aviation,  Inc. 
Los  Angeles  International  Airport 
ATTN:   Mr.  Dave  Mason 

Engineering  Data  Section 
Los  Angeles  45,  California 
M/F  Contract  AF33(038)18319 

Northrop  Aircraft,  Incorporated 
ATTN:   Northrop  Library 

Dept.  2135 
Hawthorne,  California 
M/F  Contract  AF33(600<)-22313 

Radioplane  Company 

Van  Nuys ,  California 

M/F  Contract  AF33(600)-23893 

Lockheed  Aircraft  Corporation 
ATTN:   C.L.  Johnson 
P.O.  Box  55 
Burbank,  California 
M/F  N0a(s)-52-763 

Raytheon  Manufacturing  Company 

ATTN:  Robert  Borts 

Wayland  Laboratory,  Wayland,  Mass. 

Republic  Aviation  Corporation 

ATTN:   Engineering  Library 

Farmingdale 

Long  Island,  New  York 

M/F  Contract  AF33(038)-14810 

Sperry  Gyroscope  Company 

ATTN:   Mr,  B.  Berkowtiz 

Great  Neck 

Long  Island,  New  York 

M/F  Contract  AF33(038)-14524 


Temco  Aircraft  Corp, 
ATTN:   Mr.  George  Cramer 
Garland,  Texas 
Contract  AF33(600)21714 

Farnsworth  Electronics  Co, 

ATTN:   George  Giffin 

Ft.  Wayne,  Indiana 

Marked:  For  Con;  AF33( 600) -25523 

North  American  Aviation,  Inc. 
4300  E.  Fifth  Ave. 
Columbus ,  Ohio 
ATTN:   Mr.  James  D.  Leonard 
Contract  N0a(s)  54-323 

Westinghouse  Electric  Corporation 

Air  Arm  Division 

ATTN:   Mr.  P.D.  Newhouser 

Development  Engineering 
Friendship  Airport,  Maryland 
Contract  AF33(600)-27852 

Ohio  State  Univ.  Research  Foundation 
ATTN:   Dr.  T.C.  Tice 

310  Administration  Bldg 

Ohio  State  University 
Columbus  10,  Ohio 
M/F  Contract  AF33(616)-3353 

Air  Force  Development  Field 

Representative 

ATTN:   Capt.  Carl  B,  Ausfahl,  Code  1010 

Naval  Research  Laboratory 

Washington  25,  D.C. 

Chief  of  Naval  Research 
Department  of  the  Navy 
ATTN:   Mr,  Harry  Harrison 

Code  427,  Room  2604 

Bldg,  T-3 
Washington  25,  D.C. 

Beech  Aircraft  Corporation 
ATTN:   Chief  Engineer 
6600  E.  Central  Avenue 
Wichita  1,  Kansas 
M/F  Contract  AF33(600)-20910 
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Land-Air,  Incorporated 
Cheyenne  Division 
ATTN:   Mr.  R.J.  Klessig 

Chief  Engineer 
Cheyenne,  Wyoming 
M/F  Contract  AF33(600)-22964 

Director,  National  Security  Agency 
RADE  1GM,  ATTN:   Lt .  Manning 
Washington  25,  D.C. 

Melpar,  Inc. 
3000  Arlington  Blvd. 
Falls  Church,  Virginia 
ATTN:   K.S.  Kelleher 

Naval  Air  Missile  Test  Center 
Point  Mugu,  California 
ATTN:   Antenna  Section 

Fairchild  Engine  &  Airplane  Corp. 
Fairchild  Airplane  Division 
ATTN:   L,  Fahnestock 
Hagerstown,  Maryland 
M/F  Contract  AF33(038)-18499 

Federal  Telecommunications  Lab. 

ATTN:   Mr.  A.  Kandoian 

500  Washington  Avenue 

Nutley  10,  New  Jersey 

M/F  Contract  AF33(616)-3071 

Ryan  Aeronautical  Company 

Lindbergh  Drive 

San  Diego  12,  California 

M/F  Contract  W-33(038)-ac-21370 

Republic  Aviation  Corporation 
ATTN:   Mr.  Thatcher 
Hicksville,  Long  Island,  New  York 
M/F  Contract  AF18( 600) -1602 

General  Electric  Co. 

French  Road 

Utica,  New  York 

ATTN:   Mr.  Grimm,  LMEED 

M/F  Contract  AF33(600)-30632 


D.E.  Royal 

The  Ramo-Wooldridge  Corp 

Communications  Division 

P.O.  Box  45444 

Airport  Station 

Los  Angeles  45,  California 

Motorola,  Inc. 

Defense  Systems  Lab. 

ATTN:   Mr.  A.W.  Boekelheide 

3102  N.  56th  Street 

Phoenix ,  Arizona 

M/F  Contract  NOa(s)-53-492-c 

Stanford  Research  Institute 

Southern  California  Laboratories 

ATTN:   Document  Librarian 

820  Mission  Street 

South  Pasadena,  California 

Contract  AF19(604)-1296 

Prof.  J.R.  Whinnery 
Dept .  of  Electrical  Engineering 
University  of  California 
Berkeley,  California 

Professor  Morris  Kline 
Mathematics  Research  Group 
New  York  University 
45  Astor  Place 
New  York,  N.Y. 

Prof.  A. A.  Oliner 
Microwave  Research  Institute 
Polytechnic  Institute  of  Brooklyn 
55  Johnson  Street  -  Third  Floor 
Brooklyn,  New  York 

Dr.  C.H.  Papas 

Dept.  of  Electrical  Engineering 
California  Institute  of  Technology 
Pasadena,  California 

Electronics  Research  Laboratory 
Stanford  University 
Stanford,  California 
ATTN;   Dr.  F.E.  Terman 
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Radio  Corporation  of  America 
R.C.A.  Laboratories  Division 
Princeton,  New  Jersey 

Electrical  Engineering  Res.  Lab. 
University  of  Texas 
Box  8026,  University  Station 
Austin,  Texas 

Dr.  Robert  Hansen 
8356  Chase  Avenue 
Los  Angeles  45,  California 

Technical  Library 

Bell  Telephone  Laborat-drxes 

463  West  St. 

New  York  14,  N,Y. 

Dr .  R . E .  Beam 
Microwave  Laboratory 
Northwestern  University 
Evanston,  Illinois 

Department  of  Electrical  Engineering 
Cornell  University 
Ithaca,  New  York 
ATTN:   Dr.  H.G.  Booker 

Applied  Physics  Laboratory 
Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 

Exchange  and  Gift  Division 
The  Library  of  Congress 
Washington  25,  D.C. 

Ennis  Kuhlman 

c/o  Mc  Donnell  Aircraft 

P.O.  Box  516 

Lambert  Municipal  Airport 

St.  Louis  21,  Mo. 

Physical  Science  Lab. 
New  Mexico  College  of  A  and  MA 
State  College,  New  Mexico 
ATTN:   R.  Dressel 


Technical  Reports  Collection 

303  A,  Pierce  Hall 

Harvard  University 

Cambridge  38,  Mass. 

ATTN:   Mrs.  E.L.  Hufschmidt,  Librarian 

Dr.  R.H.  DuHamel 
Collins  Radio  Company 
Cedar  Rapids,  Iowa 

Dr .  R . F .  Hyneman 

5116  Marburn  Avenue 

Los  Angeles  45,  California 

Director 

Air  University  Library 
ATTN:   AUL-8489 
Maxwell  AFB,  Alabama 

Stanford  Research  Institute 

Documents  Center 

Menlo  Park,  California 

ATTN:   Mary  Lou  Fields,  Acquisitions 


